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FREQUENCY-RESPONSE IDENTIFICATION OF XV-15 TILT-ROTOR AIRCRAFT DYNAMICS 

Mark B. Tischler, Ph.D. 

Stanford University, 1987 


The timely design and development of the next generation of tilt- 
rotor aircraft (JVX) depend heavily on the in-depth understanding of 
existing XV-15 dynamics and the availability of fully validated simula- 
tion models. Previous studies have considered aircraft and simulation 
trim characteristics, but analyses of basic flight vehicle dynamics have 
been limited to qualitative pilot evaluations. The present study has 
the following objectives: 

1. Documentation and evaluation of XV-15 bare-airframe dynamics 

2. Comparison of aircraft and simulation responses 

3. Development of a validated transfer-function description of the 
XV-15 needed for future studies 

A nonparametric frequency-response approach is used which does not 
depend on assumed model order or structure. Transfer-function represen- 
tations are subsequently derived which fit the frequency responses in 
the bandwidth of greatest concern for piloted handling-qualities and 
control-system applications. 

This study involved the planning and execution of flight tests on 
the XV-15 aircraft and piloted-simulation for four flight conditions 
from hover to cruise. Improved test techniques and pilot-training 
procedures were devised. Analytical software tools were developed (or 
adapted) which allow the identification of high-resolution spectral 
responses and the derivation and validation of multi-input/multi-output 
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transfer-function models. These techniques were applied in an extensive 
evaluation of the open-loop flight dynamics of the XV- 15 aircraft and 
simulation mathematical models. Deficiencies in the mathematical models 
were exposed and documented. Finally, a new, fully validated transfer- 
function model was derived for the hover and cruise flight conditions. 
The methods developed in this study have subsequently been applied in a 
number of other flight-test programs and have been included in the U.S. 
Army's updated helicopter handling-qualities specification. 

Approved for publication: 



Dean of Graduate Studies 
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Chapter 1 
INTRODUCTION 

1 . 1 Background and Objectives 

The tilt-rotor concept combines the hovering advantages of the 
helicopter with the cruise advantages of a fixed-wing aircraft. 
Rotor/engine nacelles at the wing tips are rotated to the vertical 
position for hovering flight and to the horizontal position for cruising 
flight. The XV-15 research aircraft (Fig. 1.1) was jointly developed by 
the U.S. Army, NASA, and the Navy to demonstrate tilt-rotor technol- 
ogy. A key objective of this project was to achieve good piloted 
handling-qualities characteristics in hovering flight by using an 
advanced stability and control augmentation system (SCAS). This objec- 
tive was emphasized because of serious handling-qualities deficiencies 
in hovering flight which were encountered with the original tilt-rotor 
demonstrator — the XV-3 (Ref. 1). Two XV-15 aircraft were developed 
under contract to Bell Helicopter Company and delivered to Ames Research 
Center (ARC) in 1980. One aircraft (N703) was retained at Ames for 
research and development testing; the other (N702) was leased back to 
the contractor for operational testing. 

Comprehensive real-time and nonreal-time simulation codes (Refs. 2, 
3) were developed to support the design and testing of the XV-15. 
Moving-base simulation facilities at ARC were extensively used for pilot 
training before the first flight tests, and subsequently for advanced 
automatic flight control system (AFCS) development. The XV-15 simula- 
tion code covers the entire operating envelope, with a full nonlinear 
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Fig. 1.1. The XV- 15 Tilt-Rotor Aircraft, (a) Hover Configuration; 
(b) Cruise Configuration. 


representation of the aircraft. Wing/body/taii aerodynamics are deter- 
mined from extensive look-up tables of full-scale wind-tunnel data 
obtained in the NASA 40- by 80-Foot Wind Tunnel facility. Rotor cal- 
culations assume quasi-steady flapping and are based on modified Bailey 
equations with uniform rotor inflow. The aerodynamic interactions 
between the two rotors, and the rotor interference with the other 
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aircraft elements are modeled in detail. Also modeled are numerous sub- 
system dynamics such as the engine drive train and governor. The XV- 15 
mathematical model was the most complex ever developed to that time for 
real-time piloted simulation at ARC (Ref. 4). The nonreal-time version 
is routinely used to support control-system development and flight-test 
planning. The XV- 15 simulation mathematical models have been exten- 
sively correlated with static trim and performance flight data; the 
comparison is generally excellent (Ref. 5). However, dynamic checks 
have been very limited (Ref. 4), with most of the validation in this 
area centered on pilot subjective comparison of the aircraft and 
motion-based simulator response. 

The author, as a staff member of an Army research team responsible 
for simulation technology, started in 1983 to conduct a comprehensive 
study to validate the open-loop dynamic response fidelity of the 
piloted-simulation mathematical model. An in-depth understanding of 
XV- 15 dynamics and the availability of fully validated simulation models 
were considered important for the timely design and development of the 
Joint Services Operational Tilt-Rotor Aircraft — the JVX, now designated 
the V-22. To fulfill these needs, the study was initiated with the 
following three major objectives: 

1. Document the open-loop dynamic characteristics of the XV-15 
aircraft from flight tests for several operating conditions including 
hover 

2. Compare aircraft and simulation response characteristics to 
identify problem areas in the mathematical modeling 

3. Develop a validated transfer-function model description of the 
XV-15 needed for future studies 
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Emphasis was initially placed on the hover flight condition, where 
unstable open-loop dynamics lead to the most critical handling-qualities 
problems . 

A key consideration in planning this study was the selection of an 
appropriate dynamics identification method. 

1.2 Dynamics Identification Methods 

Dynamics identification methods generally fall into two cate- 
gories: frequency domain and time domain. Each approach has its inher- 

ent strengths and weaknesses which make it best suited for particular 
applications. 

In time-domain (maximum-likelihood) identification (Fig. 1.2), the 
aircraft dynamics are modeled by a set of differential equations 



Fig. 1.2. Time-Domain Identification Method. 


4 










describing the external forces and moments in terms of state and control 
variables. The unknown coefficients in the equations are the stability 
derivatives, which are identified by least-squares fitting of the mea- 
sured time-responses (output-error method). Such an approach allows a 
direct comparison of stability derivatives obtained in the wind tunnel 
and those of the actual flight vehicle. Transfer functions and fre- 
quency responses may be calculated from the state-space model. A key 
aspect of time-domain identification is that an a priori model formula- 
tion must be assumed. This important step involves consideration of 
model structure, order, and important nonlinearities. Such information 
is generally not well known on a new vehicle such as the XV- 15, and 
incorrect model formulation can bias the parameter estimates (Ref. 6). 
Also, models which provide a good fit in the time-domain do not neces- 
sarily yield accurate transfer functions, since time-domain identifica- 
tion techniques weight their results more heavily at low frequency where 
most of the data points are concentrated. 

The frequency-domain identification approach shown in Fig. 1.3 uses 
spectral analysis methods to extract the frequency responses between 
selected input and output pairs. The identification results are usually 
presented in Bode-plot format, that is, log -magnitude and phase of the 
input-to-output versus log-frequency. These identification results are 
nonparametric because no model structure is assumed. As such, they are 
useful for flight-control system design and pilot-in-the-loop handling- 
qualities studies. Frequency responses obtained from real-time and 
nonreal-time simulations can be compared directly with the flight data 
to expose limitations and discrepancies in the simulator models. The 
fact that this comparison can be made initally without an a priori 
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Fig. 1 . 3 . Frequency-Response Identification Method. 


assumption of model structure or order is especially important for 
verifying mathematical models of new aircraft configurations. When the 
model structure and parametric values are required, they may be obtained 
by fitting the frequency-responses with transfer-function models to 
extract modal characteristics. Examples of this application are the 
testing of handling-qualities specifications given in lower-order equiv- 
alent system terms, and the examination of transfer function-based 
control system designs. Since this fitting procedure is completed after 
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the frequency response is extracted, the order of the transfer function 
can be selected to avoid an overparameterized model. Multi-input/multi- 
output (MIMO) frequency-response methods are suitable for extracting a 
transfer matrix which includes the important coupling effects. Finally, 
the extracted models are driven with the flight data to verify the 
time-domain characteristics. Models identified by frequency -domain 
techniques are often most accurate at mid- and high-frequency (initial 
time-history transients), which is the region of greatest concern to the 
pilot. The low-frequency and steady-state response prediction of the 
extracted models is generally not as good as in the time-domain identi- 
fication approach. 

Since the completion of the preceding objectives depends on obtain- 
ing an accurate characterization of the input-to-output transient dynam- 
ics and piloted handling-qualities of a new aircraft configuration, 
rather than on obtaining a stability derivative model (necessary for 
example to validate the wind-tunnel data base), the frequency-domain 
approach is the natural choice. 

1.3 A Historical Summary of Research in Frequency-Response 
Identification 

The earliest reported research in frequency-response identification 
of aircraft dynamics from flight-test data was conducted at the Cornell 
Aeronautical Laboratory beginning in 1945 (summarized in Ref. 7). 
Steady-state sine-wave inputs were used to (laboriously) extract the 
frequency responses of the North American B-25J (fixed-wing) aircraft. 
Then, lower-order transfer-function models were derived from a least- 
squares fit of the frequency-responses (displayed on a polar plot). 
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Fourier transform methods were subsequently developed (Refs. 7, 8) to 
allow frequency-response identification from (shorter-duration) 
discrete-maneuver data, such as that obtained from step and pulse 
inputs. These techniques were applied in flight research activities at 
the Air Force Flight Test Center (Edwards Air Force Base) during the 
1950s (see Ref. 9 for a list of references). As pointed out in Ref. 9, 
all of these early efforts in frequency-response identification suffered 
from the lack of large-scale computing power. The development of the 
fast Fourier transform (FFT) algorithms in the 1960s, and the signifi- 
cantly improved computing capabilities of this period led to much 
greater interest and success in frequency-response identification. 

A comprehensive facility for multivariable frequency-response 
(matrix) identification and analysis (FRA) based on the FFT was devel- 
oped by Twisdale (1975) at the Edwards Test Center (complete documenta- 
tion in Ref. 10). One key feature in this identification approach was 
incorporation of the ordinary, partial, and multiple coherence function 
calculations which provide an important measure of spectral estimation 
accuracy (Chap. 3). Identification was achieved from flight data of 
tracking and refueling handling-qualities tests instead of from data 
obtained with prescribed inputs. 

Marchand and Koehler (Ref. 11), of the Institute for Flight Mechan- 
ics (DFVLR) in the Federal Republic of Germany, developed a method for 
extracting the stability-derivative matrix from identified frequency- 
responses; flight data were obtained with prescribed "multi-step" 
inputs — an outgrowth of research in optimal input design (see Ref. 6). 
Frequency-response identification from flight data obtained with the 
prescribed "frequency-sweep" input was pioneered by Systems Technology, 
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Inc., (Refs. 12, 13). Advancements in lower-order transfer-function 
modeling were made by Hodgkinson et al. (Ref. 14), Bischoff and Palmer 
(Ref. 15), and Mitchell and Hoh (Ref. 16) in support of the development 
of an updated handling-qualities specification for military fixed-wing 
aircraft (Ref. 17). 

The author was the first to extensively identify frequency- 
responses and transfer-functions of rotorcraft from flight tests using 
the frequency-sweep input (Ref. 18). The frequency-domain approach used 
in this study and depicted in Fig. 1.3 draws heavily on the preceding 
researchers' efforts. The frequency-response identification (FRESPID) 
program developed for this effort (Chap. 3) was patterned after the FRA 
facility (Ref. 10). A major advancement in the FRESPID program is the 
incorporation of the chirp z-transform, an algorithm for obtaining 
frequency-responses which are of much higher quality than those obtain- 
able from standard FFT procedures. Also, the author has stressed the 
importance of deriving physically consistent transfer-function models 
and verifying those models with time-domain data not used in the identi- 
fication procedure. 

1.4 Scope of Research 

This study involved the planning and execution of flight tests on 
the XV- 15 aircraft and piloted-simulation for four flight conditions 
from hover to cruise. Improved test techniques and associated pilot 
training procedures were devised. Analytical software tools were devel- 
oped (or adapted) which allow the identification of high-resolution 
spectral responses and the derivation and validation of MIMO transfer- 
function models. These techniques were applied in an extensive 
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evaluation of the open-loop flight dynamics of the XV- 15 and simulation 
mathematical models. Mathematical model deficiencies were clearly 
exposed and documented. Finally, a new, fully validated transfer- 
function model was derived for the hover and cruise flight conditions. 

1 .5 Extensions and Applications 

As in most research efforts, the scope of this project grew well 
beyond the original objectives outlined in Sec. 1.1. The advantages of 
the frequency -domain approach for documenting the response characteris- 
tics of new configurations became readily apparent in the XV- 15 study 
and led to the use of this technique in a variety of related helicopter 
flight test projects. 

Frequency-domain testing of the Bell 214-ST single-rotor helicopter 
was completed by the author in October 1985 to support the Army's devel- 
opment of an updated handling-qualities specification (Ref. 19). A 
discussion of the testing procedures and identification results for the 
Bell 214-ST is presented in Ref. 20. This project demonstrated the 
feasibility of frequency-sweep testing and frequency-domain identifica- 
tion as a specification-compliance documentation tool. Further, the 
results showed that very low order transfer-function models can accu- 
rately predict the large motion time-domain behavior of single-rotor 
helicopters — even in hover. 

The frequency-domain identification method was also used to docu- 
ment the dynamic characteristics of the CH-47B research aircraft 
(Refs. 6, 21). Current research by the author and Kaletka (Ref. 22) 
under a memorandum-of-under standing between the United States and the 
Federal Republic of Germany is concerned with comparing frequency and 
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time-domain identification results for the open-loop XV-15 dynamics in 
hover. Research by the author and Acree (Ref. 23) concerns identifica- 
tion of XV-15 structural dynamics using the frequency-domain approach. 


1.6 Organization of Dissertation 

The central focus of this dissertation is the documentation and 
analysis of the dynamics of a unique vehicle— the XV-15. Chapter 2 
describes the general vehicle configuration and important subsystems, 
and reviews the flight test conditions. The frequency-domain method, 
briefly overviewed in the present chapter, is discussed in detail in 
Chapter 3. Chapters 4 and 5 present the dynamics identification results 
for the hover and cruise flight conditions. Frequency responses aid 
transfer-function models extracted from flight tests are compared with 
real-time and nonreal-time simulation results. Time-history verifica- 
tion responses show the ability of the extracted lower-order transfer- 
function models to accurately predict large-amplitude response charac- 
teristics of the XV-15 in hover and cruise. Chapter 6 reviews the key 
conclusions and contributions of this work. Recommendations for future 
research in frequency-response identification are discussed in Chap- 
ter 7. Appendix A reviews the important advantages of the chirp 
z-transform, as compared to the standard fast Fourier transform, for 
frequency-response identification from flight data. Appendix B dis- 
cusses the limitations in identifying open-loop vehicle dynamics from 
closed-loop flight tests. An appreciation for these limitations is 
important in applying the present identification method to future flight 
tests. The text of this dissertation draws heavily on the author's 
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publications concerning the tilt-rotor study (Refs. 18, 22, 24) and 
related research efforts (Refs. 6, 20). 


12 


Chapter 2 

XV- 15 RESEARCH AIRCRAFT AND FLIGHT TEST CONDITIONS 


An understanding of XV-15 aircraft dynamics requires an apprecia- 
tion of this unique vehicle's configuration and key subsystems. This 
Chapter first reviews the vehicle layout, operational flight envelope, 
cockpit controls, automatic flight-control system, and research instru- 
mentation, and then outlines the flight-test conditions and test inputs 
for this study. 

2.1 Aircraft Description 

2.1.1 General Layout 

The XV-15 aircraft is a lateral- tandem rotor vehicle. The rotors 
are 25 ft in diameter and the spinner-to-spinner span is 32 ft 
(Fig. 2.1). The aircraft is powered by two Lycoming T-53 turboshaft 
engines (1250 SHP each), one mounted in each wing-tip nacelle. These 
nacelles rotate from i^ = 0 deg (cruise configuration) to i^ = 90 deg 
(hover configuration). Rearward flight acceleration is enhanced by 
rotating the nacelles to i^ = 95 deg. 

The rotor system is a three-bladed prop-rotor with a stiff in-plane 
gimbal mounting to the hub. Rotor tip-path-plane orientation is con- 
trolled through standard cyclic and collective feathering of the indi- 
vidual blades. The resultant hub moments cause the entire rotor system 
to rotate (flap) as a unit, rather than each blade independently as in 
articulated helicopters. Cross-shafting between the nacelles synchro- 
nizes the rotors and provides a one-engine-out capability. An engine 


13 







governor system adjusts the pilot's power lever commands to the collec- 
tive pitch of the blades to maintain constant rotor rpm. 

2.1.2 Operating Flight Envelope 

The operational envelope shown in Fig. 2.2 is from slow rearward 
flight to a maximum true airspeed of V t = 300 knots (with the nacelles 


NORMAL RATED 
POWER LIMIT 



Fig. 2.2. Airspeed-Altitude Envelope. 

locked at i^ = 0 deg). When the nacelles are not in their locked 
position, the maximum airspeed is restricted to V fc = 170 knots because 
of aeroelastic stability limitations. The "conversion corridor" of 
Fig. 2.3 defines the allowable range of nacelle-angle/airspeed 
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0 FLIGHT TEST CONDITIONS 
FOR PRESENT STUDY 



Fig. 2.3. Conversion Corridor. 

combinations in the transition flight regime between hover and cruise. 
The lower boundary on nacelle angle is determined by wing-loading 
limits, and the upper boundary is determined by blade loading limits. 

The nominal weight of the XV- 15 is 13,000 lb and the flight endurance is 
approximately 1 hr. 

2.1.3 Vehicle Control 

The pilot's cockpit controls in the XV- 15 are the same as those 
found in conventional helicopters. A center stick is used to control 
the pitch and roll motions. Pedals are used to control yaw motions, and 
a power lever is used to control vertical motion in hover and airspeed 
in forward flight. The means with which these cockpit controls generate 
the required forces and moments depends on the aircraft configuration. 
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The control method for helicopter-mode flight is shown in 
Fig. 2.4. Rolling moments are generated with differential collective, 
pitching moments with uniform longitudinal cyclic (fore-aft tip-path- 
plane rotation), yawing moments with differential longitudinal cyclic, 
and heave forces with uniform collective. In cruising flight, the rotor 
controls are phased out and control moments are generated with standard 
aircraft aerodynamic surfaces; the ailerons produce rolling moments, the 
elevator produces pitching moments, and the rudder produces yawing 
moments. 





Fig. 2.4. Control Method for Helicopter-Mode Flight. 


Speed control is achieved through the throttle/governor system. 

The rotor control (swashplate) and aerodynamic surfaces receive commands 
from the cockpit via mechanical linkages and hydraulic actuators. These 
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commands are augmented by signals from the automatic flight-control 
system through series actuators. 

The phase-out of the rotor controls from hover to cruise is auto- 
matically scheduled as a function of nacelle angle i N . The aerodynamic 
surfaces are actuated throughout the entire flight envelope, although 
they have no effect at very low airspeeds. It is most expedient to 
refer all of the open-loop vehicle response characteristics (e.g., 
frequency responses, transfer functions, etc.) to these surface deflec- 
tions since, neglecting the small servo lags, these are related to the 
sum of the pilot and SCAS inputs (total input to the aircraft) through a 
mixing ratio which is constant across the entire flight envelope. 

In the helicopter configuration, the XV-15 open-loop dynamics are 
typical of hovering vehicles, although the number of cross-coupling 
paths is much lower than in single-rotor helicopters. The planar sym- 
metric, lateral-tandem configuration yields vehicle dynamics which are 
essentially decoupled between the longitudinal and lateral/directional 
degrees-of-freedom. One-way coupling from roll input to yaw response is 
significant. This is due to the differential rotor torque which accom- 
panies the differential collective inputs for roll control. The open- 
loop pitch and roll dynamics in hover are characterized by highly 
unstable low-frequency attitude-speed divergences. The heave and yaw 
dynamics are decoupled from the attitude motions and are essentially 
first-order. The time-constants of these dynamics are very long 
(= 10 sec) because the tilt-rotor configuration has a relatively high 
disk-loading (low heave damping) and no tail-rotor (low yaw damping). 

Limited evaluations of the SCAS-off handling-qualities in hover 
have been conducted (Ref. 25). These evaluations indicate that the 
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severe attitude instability and small translational damping result in 
very poor handling-qualities. Precision hover tasks (out of ground- 
effect) were rated with Level II handling-qualities ("adequate perfor- 
mance achieved with high level of pilot workload"). Precision transla- 
tion tasks were rated with Level III handling-qualities ("unsatisfactory 
performance"). Quantitative documentation of the open-loop XV-15 dynam- 
ics in hover is critical to "convert" these pilot-opinion ratings into 
engineering requirements for future tilt-rotor configurations. 

2.1.4 Automatic Flight Control System 

A nominal stability and control augmentation system (SCAS) was 
incorporated in the XV-15 to enhance the poor bare-airframe dynamic 
characteristics. Improved closed-loop handling qualities were achieved 
with an advanced SCAS developed by Churchill and Gerdes (Ref. 25) 
between 1982 and 1984. This advanced system combines feed-forward 
shaping, model-following, and disturbance rejection to achieve crisp 
first-order responses in the attitude rates and improved inherent sta- 
bility. The closed-loop response time constants are about 0.4 sec in 
pitch and roll and about 0.6 sec in yaw. An attitude-retention system 
based on attitude-angle feedback provides hands-off stability when the 
controls are in their neutral position. In forward flight, the control 
system provides constant load factor responses to longitudinal stick 
inputs and coordinated turn responses to lateral stick inputs. These 
improved characteristics give the aircraft Level I handling qualities 
("performance satisfactory without improvement") for the precision 
piloting tasks (Ref. 25). 

The SCAS is a two-channel (summed) fail/operate system. Thus, 
failure of either channel results in reversion to a single (lower-gain) 
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channel. Failure of the second channel reverts the aircraft back to the 


SCAS-off configuration. The SCAS authority is limited to 20 % of maximum 
stick deflection in the pitch axis, and 30 % of the maximum control in 
roll and yaw. 

2.1.5 Research Instrumentation 

The XV- 15 is heavily instrumented to provide real-time and post- 
flight engineering data. Measurements of over 150 variables including 
angular rates, attitudes, accelerations, surface deflections, cockpit 
controls, and structural loads are recorded on board at 250 Hz on 
pulse-code-modulation (PCM) tapes. Besides this on-board recorder, a 
telemetry (TM) link provides safety-related data which are monitored by 
a six-member flight-test ground crew. This extensive instrumentation 
system is carefully maintained and regularly calibrated. After the 
flight tests, the PCM tapes are digitized, and the data are converted to 
engineering units. 

The three-axis gyroscope package has a bandwidth of about 20 Hz, 
which is well beyond the frequency range of interest for rigid-body 
dynamics identification. Analog measurements of controls and responses 
are conditioned with 50 Hz analog pre-filters to reduce digital alias- 
ing, while maintaining a broad band of accurate dynamic measurements. 
Matching the filters on the input and output signals minimizes the phase 
distortions in the identified frequency responses. 

2.2 Flight Test Conditions 

The following four flight test conditions were selected to span the 
full operating range of the XV- 15 aircraft: 
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1. Hover: Ambient wind less than 5 knots, i N = 90 deg, gear down, 

out-of-ground effect (altitude = 100 ft) 

2. Low-speed transition: Indicated airspeed V^ = 100 knots 

(calibrated airspeed, V c = 105 knots), i N = 70 deg 

3. High-speed transition: Indicated airspeed = 120 knots 

(calibrated airspeed, V c = 128 knots), ij| = 30 deg 

4. Cruise: Indicated airspeed = 170 knots (calibrated air- 

speed, V c = 180 knots), i N = 0 deg 

These flight conditions are noted on the conversion corridor of 

Fig. 2.3. 

Dynamics identification tests were conducted from February through 
December 1983, in about 12 flight hours. During the same period, 
frequency -sweep tests were conducted in the Vertical Motion Simulator 
(VMS) to document the real-time XV- 15 mathematical model. Subsequent 
analyses of aircraft and simulation flight dynamics concentrated on the 
hover and V^ = 170 knots cruise conditions. 

Frequency-sweep (Sec. 3.1.2) and step inputs were executed in each 
axis for all four flight conditions. In hover, three frequency-sweeps 
were conducted to ensure that sufficient dynamic data were obtained for 
good identification of this most important and (difficult to analyze) 
condition. In the remaining flight conditions, where the vehicle is 
much more stable, only two frequency-sweeps were required. 

The high degree of open-loop pitch and roll instability in the 
hover flight condition dictated that longitudinal and lateral stick 
frequency-sweeps in this flight condition be conducted with all SCAS 
channels ENGAGED. Pedal-sweeps were conducted with the yaw-SCAS disen- 
gaged, because yaw-SCAS failures occurred in mid-run resulting from the 
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relatively large angular rates which were generated. The decoupled and 
stable nature of the yaw (and heave) dynamics allowed SCAS-off sweeps to 
be conducted in these axes without difficulty. Step inputs were exe- 
cuted in the SCAS-on and SCAS-off configurations. SCAS-off step 
responses were completed to verify that the extracted transfer-functions 
reflect the open-loop response characteristics and not those of the 
inverse feedback dynamics (Chap. 3). These inputs were completed with 
the off-axis channels engaged (e.g., open-loop pitch-axis inputs were 
conducted with roll and yaw channels engaged), to reduce the level of 
coupled, unstable dynamic response. SCAS-on step responses were also 
completed because the SCAS helps to steady the initial conditions, which 
is helpful in exposing small differences between the model and aircraft 
dynamics. 

In the transition and cruise flight conditions, the longitudinal 
frequency-sweeps were conducted with pitch SCAS engaged . Frequency- 
sweeps in the lateral and directional axes were completed with the roll 
and yaw SCAS channels disengaged because initial SCAS-on results exhib- 
ited unacceptable levels of aileron/rudder correlation (Chap. 3). Step 
inputs in all axes were completed for both SCAS-on and SCAS-off config- 
urations as in hover. 

The next chapter discusses the dynamics-identification method in 
detail. 
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Chapter 3 

FREQUENCY-DOMAIN IDENTIFICATION METHOD 

The frequency-domain method of Fig. 1.3 was briefly outlined in 
Chapter 1. In this Chapter, the individual analysis steps of the 
approach are discussed in detail. 

3.1 Flight Test Technique 

3.1.1 General Requirements 

The principal ingredient for success in any identification scheme 
is the selection of an input signal that excites the vehicle in all of 
its dominant modes of motion in the frequency range of interest. In 
piloted handling-qualities studies, the important frequency range is 
approximately 0. 2-6.0 rad/sec (Refs. 17, 19). To achieve good 
frequency-domain identification in this range, the magnitude of the 
input autospectrum should be nearly constant. From a piloting stand- 
point, the selected inputs should be easily repeatable and not involve 
drastic maneuvers or significant changes in the equilibrium flight 
condition. One input which fulfills the preceding requirements and 
which has been used successfully in a number of rotorcraft flight tests 
(Refs. 6, 18, 20, 21, 26) and non-rotorcraft flight tests (Refs. 12, 13, 
27) is the "frequency-sweep." 

3.1.2 Frequency Sweep 

A typical lateral-stick (6 LflT ) frequency-sweep completed during the 
XV-15 hover tests and developed by the author is shown in Fig. 3.1. 
Pilot-generated inputs were used instead of computer-generated inputs. 
The sweep is initiated with two low-frequency sinusoid-shaped cycles, 
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FIRST SECOND 



with periods corresponding to the lower bound of the desired identifica- 
tion frequency range. These cycles ensure good excitation of the low- 
frequency dynamics. The desired lower bound of u> min = 0.2 rad/sec 
gives a low-frequency period of about T max = 30 sec. However, attempts 
to execute this low-frequency input in flight resulted in undesirably 
large vehicle motions. Therefore, it was necessary to reduce the period 
of the low-frequency cycle to T max = 20 sec (Fig. 3.1). The magnitude 
of the control input is adjusted to keep roll and pitch attitudes, as 
well as resulting translations, comfortable and reasonably close to 
their trim values. The lower pitch and roll damping in hovering flight 
demands smaller pilot inputs than those which are used in transition and 
forward flight. 

After the initial two low-frequency cycles, the control is moved at 
gradually increasing frequency for another 50 sec. This increasing- 
frequency (as opposed to a decreasing-frequency) input allows the tran- 
sients of the low-frequency modes to persist for a few more cycles 
during the remainder of the run, thus improving the low-frequency 
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identification; the high-frequency modes need less excitation time to 
persist the same number of cycles, so they are excited toward the end of 
the input. By the end of the run, the control is being driven quite 
rapidly (about 4 Hz in Fig. 3.1) with generally smaller displacements. 
The control is then returned to trim, ending the approximately 90-sec 
test period. The selection of the (rather long) overall 90-sec test 
period was found (empirically) to 

1 . Allow a good identification of the low-frequency modes 

2. Result in an even pilot-excitation of the vehicle dynamics over 
the identification frequency range 

3. Yield a large quantity of data, thereby reducing the variance 
of the spectral estimate (see Sec. 3.2.3) 

At least two sweeps are performed consecutively for each control to 
ensure a sufficient amount of good quality dynamic data for the spectral 
analysis. Frequency-sweeps are conducted individually on each of the 
four primary aircraft controls (lateral and longitudinal stick, pedals, 
and power lever) for one flight condition, before other conditions are 
investigated. 

3.2 Frequency-Response Identification 

The most important step in the frequency-domain identification 
method is the accurate estimation of input-to-output spectral responses, 
since all succeeding calculations are based on these results. Consider- 
able effort was invested at the beginning of this study in evaluating 
and implementing various spectral-analysis methods. Numerical studies 
based on reconstructing frequency responses from simulated time- 
histories of known transfer functions were useful in evaluating the 
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various computational methods. However, methods which proved successful 
with simulated data did not always give satisfactory results when 
applied to flight data. Good quality results at low frequency initially 
proved very difficult to obtain. This was due to the poor spectral 
content which generally occurs at low frequency because of the lower 
pilot input amplitudes (see Fig. 3.1). Considerable improvements were 
achieved with more advanced analysis methods developed as a result of 
identification efforts on the XV- 15, CH-47, and Bell 214-ST aircraft 
(Refs. 6, 18, 20-24). These methods draw heavily on the signal process- 
ing literature (Refs. 28-32). 

Figure 3.2 shows the computational procedures which are performed 
by the FRESPID program to generate the spectral responses of selected 
input and output time-histories. This program was jointly developed by 
the author and J. G. M. Leung. t The FRESPID program operates in a batch 
processing mode, under the control of a command file. As illustrated in 
Fig. 3*3, this file contains command directives and flight-data record 
identifiers which specify the desired input/output records and data 
processing options. The next two sections describe in more detail the 
computations which are performed by the FRESPID program. 

3.2.1 Data Preparation 

As indicated in the flowchart of Fig. 3.2, a large part of the 
analysis involves data preparation before applying the Fourier transform 
algorithm. These preparation steps are crucial to obtaining an accurate 
frequency-response estimate. 


t Electronics Engineer; Flight Experiments Support Branch, NASA Ames 
Research Center, Moffett Field, Calif. 
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Fig. 3.2. Flowchart of Computational Procedures Performed by the FRESPID Program. 











Command file directives 


Remarks 


V=0. P/DAIL 883,884 60S WINDOW [OC128] 
WINDOW 

MYDISK: [ TISCHLER . TILTBFILE ] XV15W60S 
LINK 

D645R883. 150 V012R883.150 
D645R884. 150 V012R884.150 
FILTER 

MYDISK: [TISCHLER. BFILE] AALIAS .62H 
0 . 1 , 10.0 
0C128 
STOP 


command file identifier 
key word: window specification 
window specification file 
key word: file concatenation 
input file-1 output file-1 
input file-2 output file-2 
key word: digital filtering 
impulse response file 
Chirp z-transform freq. range 
output file name 
program terminator 


Fig- 3-3- Sample Command File for FRESPID. 


The dc (offset) components and linear drifts (trends) of the 
selected input and output time-histories are first removed to prevent 
oscillation in the spectral calculations. Since the frequency-sweep 
begins and ends in trim and is a roughly symmetric excitation 
(Fig. 3.1), the dc and drift are simply determined from the mean offset 
and slope of the time-history trace (Ref. 28). When (nr) multiple runs 
of good spectral quality are available, they are concatenated to form 
composite input and output time-histories. The concatenated time- 
histories are then digitally filtered to eliminate the jitter (digital 
noise) which arises from asynchronous sampling of the instruments 
("sample-skewing"), and can cause distortion in the spectral estimate. 
The filter cutoff frequency is 62.5 Hz, which is 50 % of the Nyquist 
frequency. Finite- impulse-response (FIR) filters are used [instead of 
recursive infinite- impulse-response (HR) filters] because they have 
magnitude characteristics that can be easily tailored, and they exhibit 
a linear phase shift with frequency (independent of the magnitude 
response) (Ref. 28). The additional computational time for the FIR 
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filters compared to that of IIR filters is not a concern for off-line 
analyses. 

The filtered concatenated records are partitioned into (K) tapered 
and overlapping subrecords of (L) discrete data points. These sub- 
records, or time-history sections, are often referred to as "windows," 
because the data contained therein are the only part of the time-history 
that the fast Fourier transform (FFT, Sec. 3.2.2) can "see" at one 
time. Applying a cosine-squared tapering (weighting) function (also 
referred to as "Hanning window") prevents side lobes and leakage which 
would otherwise occur with a rectangular (nontapered) window. Over- 
lapping the windows makes more efficient use of the data at the edges of 
the tapered window, thereby reducing the spectral bias and variance 
compared with the non-overlapped tapered estimate (Ref. 28). 

The width of the window (LAt, sec; At is the time increment) is an 
important parameter in determining the quality of the spectral identifi- 
cation. The windows are optimized for each run (and flight condition) 
to achieve maximum spectral accuracy over the desired frequency range. 

As discussed in Sec. 3.2.3, the random error is inversely proportional 
to the square root of the number of (independent) time-history sec- 
tions— which is the effect of averaging. So for a given total run time, 
smaller windows yield more averaging and lower spectral variance. How- 
ever, this reduction in window length also reduces the low-frequency 
content in each window and, accordingly, in the overall identifica- 
tion. Thus, the selection of window length involves a compromise 
between a high number of averages and adequate low-frequency signal 
content. By concatenating good repeat runs into a single time-history, 
longer windows can be selected without compromising the number of 
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averages. This yields wide-band identification and low spectral 
variance. When only one suitable run is available, shorter windows are 
necessary to maintain a sufficient number of averages, and the low- 
frequency identification is compromised. A window length which contains 
roughly two periods of the lowest-frequency wave of interest 
(LAt = 2T mav = 40 sec) is generally a good compromise. 

3.2.2 Spectral Functions^ 

The finite Fourier transform X(f,T) of a (finite length, T) con- 
tinuous time-history subrecord x(t) (e.g., one window of time-history 
data) is: 


X(f,T) 


f T *(t)e-J 2 " ft dt 

J o 


(3.1) 


When the time-history data are from a discrete sequence x n , the 
Fourier transform integral relation is approximated by the discrete 
Fourier transform (DFT): 


N-1 

X(f ) = X(kAf) = At £ x exp[-j2*(kn)/N] ; 

n=0 


k = 0, 1, 2,..., N-1 

(3.2) 


where 

X(f^) = Fourier coefficients 
x n = x(nAt) 

At = time increment 

N = number of discrete frequency points 


t Background information on spectral analyses and the chirp z-transform 
is given in Appendix A. 
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For the DFT, 


N = L = number of time-history points in the sequence, x n (3.3) 


so, 

T = LAt = time duration of sequence (i.e., "width of window") 


Af = 


= frequency resolution, Hz 


(3.4) 


p 

The evaluation of Eq. (3.2) requires approximately multiply-add 

operations for each window of data. 

The fast Fourier transform (FFT) refers to efficient algorithms 
(there are many, Ref. 28) for determining the DFT of Eq. (3.2). The FFT 
procedures require approximately 4N log2 N computations, which gives a 
speed ratio of (Ref. 28): 

N 2 

speed ratio = -rrr-T ^ = 160 for the present analysis 

AN log2N 

indicating a considerable computational savings. One particularly 
flexible FFT algorithm is the chirp z-transform (also referred to as 
the "zoom transform" or CZT). This algorithm (Refs. 30, 32) allows the 
rather severe restriction of Eqs. (3.3) and (3.4) to be relaxed, which 
yields much finer frequency resolution (Af) for a given window size 
(LAt). Another key advantage of the chirp z-transform is that the 
number of time-history points L can be arbitrary, whereas L must be 
a composite integer (an integer power of 2) in standard FFT proce- 
dures. These capabilities of the chirp z-transform allow the extraction 
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of high-resolution spectra in a narrow frequency band, and an increase 
in the identified dynamic range, especially at the low-frequency end 
(Ref. 32). The chirp z- transform has been found in this and in a number 
of other studies by the author to be well suited for frequency-response 
identification from flight data (Refs. 6, 20, 21, 23). 

The one-sided input autospectral-density estimate for the sub- 
record x n is determined from the chirp z-transform (Fourier) 
coefficients: 


^xx (f k ) = lu l X(f k } ! 2 * k = 0, 1,..., N/2 (3.5) 

where U is the scale factor for window tapering [U = (8/3)^^ for 
Hanning window, Ref. 28], 

The one-sided output autospectral-density estimate is similarly 
obtained from the output subrecord by replacing y (output record) for 
x in Eqs. (3. 1)-(3.5): 


5 yy (f k ) = lu l Y(f k ) ! 2 ’ k = 0, 1,..., N/2 (3.6) 


The one-sided cross-spectral density estimate is determined byt : 


6 xy (f k } = fu [ x *<V Y ( f k >] > k = 0, 1,..., N/2 (3.7) 


t 


Asterisk (*) denotes complex-conjugate. 
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Finally, the total spectral-function estimates for the entire 
(concatenated) time-history is obtained from a linear average of the 
spectra for the (K) overlapped subrecords: 


Gxx ( f k ) = avera 6 e d i n P ut autospectral-density estimate 
G (f^) = averaged output autospectral-density estimate 

J J 

A 

G X y(fj £ ) = averaged cross-spectral density estimate 


These quantities are physically interpreted as giving the 
distribution of the mean-squared response of the respective signals 
(i.e., xx, yy, xy) as a function of frequency. For illustrative 
purposes, the spectral-density magnitudes are presented in power- 
dB (= 10 log 10 |G xx (f k ) | , etc.) which gives the distribution of the root- 
mean-squared (rms) response. These results are the bases for the 
frequency-response calculations in the next section. 

3.2.3 Single-Input/Single-Output Frequency-Response Calculations 
in the Open-Loop 

Once the input, output, and cross-spectral density estimates 
[(G xx (f k )» G yy (f k ), G X y(f k ), respectively] have been determined for a 
selected time-history pair (x, y in Fig. 3. 4) the single-input/single- 
output frequency response is estimated by 


x (t) 


H (f) 


y (t) 


Fig. 3.4. Single-Input/Single-Output Open-Loop Response 
Identification. 
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component of the Fourier series). The remaining output harmonics are 
called remnants. Most of these are largely filtered out by the low-pass 
nature of the aircraft dynamics. 

The piloting technique and time-history concatenation methods 
described above result in input amplitudes that are not constant over 
the frequency range. Therefore, the resulting describing function is 
representative of the average input amplitude over the ensemble time- 
history. Since the pilot inputs are typical of those experienced during 
normal flight operations, this approximation produces a satisfactory 
description of the vehicle dynamics in its normal operating environment. 

A good indication of the quality of the first harmonic as a model 
of the input-to-output dynamics is obtainable from the coherence func- 


tion 



) defined by 
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2 


(3.9) 


which may be interpreted as that fraction of the output spectrum that 
can be accounted for by linear relation with the input spectrum 
(Ref. 29). When the process under investigation is perfectly linear and 
the spectral estimates are noise free, the coherence function will be 

unity within the excitation frequency-range. (The meaningless result 

2 

of y =1 also occurs when only one window section is used — that is, 

xy 

no averaging). 

A coherence function of less than unity results from three basic 
causes (Refs. 28, 34). First, nonlinearities in the system may produce 
remnants which are not accounted for by the first harmonic approxima- 
tion. When the input and output excursions are small, the system non- 
linearities are less important. This is generally true for the cruise 
flight conditions. However, in hover, where low-frequency inputs induce 
large motions relative to the trim condition, nonlinear effects become 
more important; a significant reduction in the coherence function will 
result when first harmonic approximations are inadequate for represent- 
ing these nonlinearities. In the higher-frequency range, where the 
excursions are smaller, the coherence function again approaches unity. 

A second and very common cause of reduced coherence is input and 
output noise. The effect of noise sources differs chiefly depending on 
whether they are correlated with the input and output signals, and 
whether they act on the "process" or on the "measurement." Process 
noise (e.g., turbulence) contributes to the aircraft response, and so 
acts the same as a secondary (but unmeasurable!) input; the problem of 
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secondary inputs is discussed below. Measurement noise does not con- 
tribute to the aircraft response, but it does contaminate the recorded 
input and output signals. Process and measurement noise are important 
because they cause errors in the estimate of the frequency response. 
These errors are categorized as "bias” (systematic) or "random" (nonsys- 
tematic). Bias error exists when the expected value of the ensemble- 

A 

averaged frequency-response estimate E [ H( f ) ] is different from the 
actual frequency response H(f) or 


b[H(f ) ] = E[H(f ) ] - H(f ) 

Random error exists when there is a dispersion of the frequency-response 
estimates about the expected value for a collection of repeated runs: 

o[H(f)] = + f^E[(H(f) - E[H]) 2 ] 

Measurement noise is generally considered to be uncorrelated with 
the measured signal. Uncorrelated measurement noise at the output [v(t) 
in Fig. 3.5] does not bias the spectral estimate of Eq. (3.8). However, 
uncorrelated measurement noise at the input [u(t) in Fig. 3.5] does bias 
the estimate as a function of the noise-to-signal ratio (Ref. 28); so, 



Fig. 3.5. Uncorrelated Measurement Noise at the Input and Output. 
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input sensors must be of the highest possible quality (ideally noise- 
free). Correlated measurement noise at the input or output or both will 
also bias the estimate. Uncorrelated measurement noise also causes 
random error, which is reflected by the drop in the coherence func- 
tion. This coherence drop is often noticeable in the hover flight 
condition, where the pilot's low-frequency inputs are generally small. 
This also often occurs at higher frequency (above about 1 Hz), where the 
aircraft response is generally small. For a given level of coherence, a 
significant improvement in the spectral accuracy is obtained when multi- 
ple runs are concatenated, since the random (normalized) error 
(e r = o[H]/H) is inversely proportional to the number of independent 
time-history averages (Ref. 34): 


e 


r 



(t) 


( 3 . 10 ) 


and 


n d 



= number of independent time-history averages 


where 

Tp = Total concatenated run length, sec 
L = Number of time-history points in the window 
At = Time increment 


t The factor of f^O.55 accounts for the 45% reduction in spectral 
variance due to the 50 % window overlapping used in this study (Ref. 31). 
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Based on the average window length 2T max = 40 sec, and two concatenated 

90-sec runs, a random error of e r < 0.1 (less than 10^) is achieved for 

2 

a coherence function of y >0.85. 

x y 

The third source of reduced coherence results from secondary 
inputs; these are defined as inputs (other than the control which is 
undergoing a frequency-sweep, x in Fig. 3.6) which excite the system. 



SECONDARY 

INPUTS 


Fig. 3.6 Effect of Secondary Inputs on Single-Input/Single-Output 
Identification. 

Secondary inputs occur from gusts and turbulence (process noise), and 
from off-axis control activity of the pilot or SCAS. The gust and 
turbulence inputs cannot be measured, they contribute to the measured 
response of the system, and they generally evoke correlated regulatory 
inputs from the pilot and the SCAS; therefore, the extracted frequency- 
responses are biased (see Sec. 3.2.4). Furthermore, the random error 
increases. Therefore, frequency-sweeps must be conducted in conditions 
of minimum wind and turbulence. The author's experience indicates that 
hover tests must be conducted with steady winds that are no greater than 
5 knots and that cruise tests must be conducted with turbulence that is 
no greater than ±1 knot. The residual effects of small levels of gusts 
and turbulence are minimized by the concatenation of multiple frequency- 
sweep records. Pilot or SCAS inputs or both for regulation against 
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excursions which occur in the coupled degrees-of-freedom also contribute 
to the system response; however, these inputs are measurable. If the 
secondary and primary inputs are not fully correlated , their effects can 
be separated by using the multi-input/multi-output identification method 
described in Sec. 3.2.5. 

3.2.4 Identification of the Open-Loop Single-Input/Single-Output 
Frequency Response From Closed-Loop Tests 
When the bare-airframe characteristics of the vehicle under study 
are highly unstable, as are those of the XV-15 in hover, execution of a 
90-sec frequency-sweep on the open-loop configuration is not practi- 
cal. The single- input/single-output bare-airframe frequency-response 
can be extracted from the closed-loop flight data using Eq. (3.8) sub- 
ject to an important qualification. Referring to Fig. 3.7, the closed- 
loop system responds to commanded inputs (^AT^ and P rocess noise (n). 
Identification of the open-loop frequency response (p/$ a ) requires the 
measurement of the output (p) and the surface deflection (S_) 

a 



(SCAS) 


Fig. 3.7. Single Degree-of-Freedom Closed-Loop Roll-Response Model. 


39 



signals. This total surface deflection is made up of components from 
the pilot (« LAT ) and fche SCAS feedback Uf )• Thus the identified 

el 

frequency response of p/6 a is also expressible as 


When, 



P 

6 LAT " pG F 

P 



>> 6 


LAT’ 



(3.11) 


( 3 . 12 ) 


which is the troublesome result of identifying the inverse feedback 
frequency-response “1/Gp^ in closed-loop pilot tracking tasks without 
a significant known external input (Ref. 35). The limiting condition of 
Eq. (3.12) is reachable only in the case of system excitation by process 
noise (n) (e.g., turbulence) with no pilot inputs ($LAT = However, 

biases in the open-loop frequency-response estimate obtained from 
Eq. (3.8) will occur even with non-zero pilot input ( fi LAT * whenever 
process noise is present; this occurs since the feedback loop causes the 
surface deflection to be partially correlated with the process noise 
(Appendix B). The amount of bias error depends on the noise-to-signal 
rms ratio o(w) at each frequency; the numerical study presented in 
Appendix B shows that the (normalized) bias errors become significant 

A 

(e b = b[H]/H > 0.10; greater than 10J) when the system rms excitation 
that is due to process noise is greater than 33% of that that is due to 
pilot's inputs [o(a>) > 0.33]. Total error (composed of normalized 
random and bias error) becomes significant (greater than 10% ) at lower 
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noise-to-signal rms ratios [3(u>) > 0.12] depending on the level of 
feedback gain. 

Therefore, the key requirement for identifying the open-loop fre- 
quency response from closed-loop flight data is that the total surface 
deflection ( 6 a ) must contain a significant component from the pilot 
( ^LAT^ which is uncorrelated with the output in the frequency range of 
interest. The "programmed" frequency-sweep input has good input power 
across the entire frequency range of concern, and persists throughout 
the duration of the time-history; as a result, it is especially good for 
closed-loop testing. Multi-step inputs are not particularly well suited 
for open-loop identification in closed-loop tests because the excitation 
occurs principally at the beginning of the run. Most of the long-term 
response is excited by fully-correlated inputs from the SCAS. The use 
of the frequency-sweep input and conducting the flight tests in low- 
turbulence conditions ensures that an accurate open-loop identification 
is achieved. SCAS-off step responses of the vehicle are compared with 
the extracted transfer-function model responses to verify that the true 
(unbiased) bare airframe frequency response is identified. 

Identification of the unstable bare-air frame dynamics is thus 
greatly simplified since the pilot is flying the closed-loop (stable) 
aircraft. When turbulence is present, the remaining (secondary) 
degrees-of- freedom will be excited. For example, pitch and yaw motions 
will generally be present during lateral stick frequency-sweeps. How- 
ever, if the roll response is completely decoupled from the pitch and 
yaw responses, the identified roll frequency-response and coherence 
functions will not be affected. For such aircraft, the pilot or 
stability-augmentation systems should be used to regulate against 
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introducing disturbances into the secondary degrees-of-freedora. A 
problem arises, though, when coupled aircraft dynamics are considered. 

3.2.5 Multi-Input/Multi-Output Frequency Response Identification 
in the Closed-Loop 

Consider the case of coupled roll-yaw dynamics with SCAS feedbacks 

in both channels, as illustrated in Fig. 3.8. Following the preceding 

discussion, we assume that the total surface deflections (6 & and 6 r ) 

contain substantial components from the pilot inputs an ^ ^ped^ 

which are uncorrelated with the output responses. 

Now suppose that a lateral-stick sweep is being executed to obtain 

the open-loop transfer function for yaw response to aileron (r/6 a ). 

Unlike the previous decoupled case, aileron inputs induce coupled 

responses in the yaw degree-of-freedom. If SCAS and pilot yaw feedbacks 

(G R ) act to restrain these coupled yaw excursions, the identified 
‘‘r 

transfer function r/6 0 will be a greatly distorted indication of the 

a 

actual bare-airframe transfer function. 

In a general sense, yaw-rate feedback decouples the yaw and roll 
dynamics^ ; this may be desirable for assessing the open-loop roll 
response (p/6 a ) using the simplified one degree-of-freedom model 
(Fig. 3.4). However, when the interaxis-coupling characteristics are 
sought (p/5 r , r/6 a ), artificial decoupling is not satisfactory. In 
hovering rotorcraft, rotor flapping dynamics induce interaxis coupling 
effects which may degrade piloted handling-qualities (Ref. 36); there- 
fore, the identification of interaxis coupling responses is a specific 
concern . 


t This assumes that is small, as it is in the present XV- 15 case. 
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Fig. 3.8. Coupled Roll/Yaw Response Model. 


The two- input/s ingle-output identification problem illustrated in 
Fig. 3.8 is a special case of the more general multi-input/single-output 
analyses discussed in Refs. 10, 28, and 29. The present analysis is 
limited to the two- input problem since the side-by-side rotor configura- 
tion of the XV-15 aircraft induces coupling between the roll and yaw 
degrees-of-freedom only (Chap. 2). (Identification of open-loop 
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dynamics for single-rotor helicopters may require consideration of more 
than two-control inputs . ) 

The components of the yaw-rate spectrum resulting from the aileron 
and rudder inputs are determined from 


V * fe) \«a * fc) °Vr 

V * fe) V. * fe) \‘ r 


(3.13) 


(3.14) 


The desired frequency responses are yaw-rate response to aileron (r/6 a ) 
and yaw-rate response to rudder (r/« r ). These are given by the simulta- 
neous solution of Eqs.(3.13) and (3.14): 


G 4 r' 1 - <G S 5 G « r /G 4 4 G 4 r )] 
a a r r r r a 


°4 4 <’ 
a a 


A . > 

a r 


(3.15) 


G 4 r' 1 - <°4 4 G « r /G 4 4 G « r )] 
r _ r r a a a a r 

6 . 2 , 

r S 5 " ^{ 5 

r r a r 


( 3 . 16 ) 


It is clear from the above equations that the aileron and rudder 
inputs must not be fully correlated since this would result in an input 

p 

cross-coherence (y. . ) of unity, thereby rendering the desired fre- 
6 6 
a r 

quency responses undefined (the "multipath propagation problem”, 

Ref. 34). Then, only the (assumed) decoupled roll-rate response to 
aileron can be identified. Note also that if the inputs are fully 
uncorrelated . Eqs. (3.15) and (3.16) reduce to the single-input/single- 
output relations from Eq. (3.8). 
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Equations (3.15) and (3.16) can be written more compactly by 


defining 


V'V 2 V ' ,G ‘ * °« r /G « . G « r» 

ar a arr rra 


(3.17) 


G = G (1 - y ) 

6 6-6 ' 6 6 v { ' 
a a r a a a r 


( 3 . 18 ) 


; 5 r-6 2 G s r [ ’ - (G « « G i r /G S « G « r » 

ra r raa aar 


(3.19) 


6 6 ■ 
r r 


sG { « <’ - 
r r 


A, 1 

a r 


( 3 . 20 ) 


Each quantity defined in Eqs. (3- 17)-(3-20) has important physical 
significance. For example, G is the cross-spectrum from aileron 

o V • o 

a r 

input to yaw-rate response (6 a to r) with the linear effect of yaw-rate 
due to rudder excitation removed. This quantity is termed the "condi- 
tioned (or residual) cross-spectrum" because the "ordinary" cross - 

spectrum (G^ ) has been conditioned to remove the biasing effects of 

°a r 

partially-correlated rudder excitation. Similarly, the quantity 

G is termed "conditioned input autospectrum," because the ordi- 

a a r 

nary input autospectrum (G r ) has been conditioned to remove the 

o o 
a a 

partially-correlated contribution of 6 r inputs. The remaining quanti- 
ties of Eqs. (3.19) and (3.20) have analogous interpretations. 

The desired frequency responses of Eqs. (3.15) and (3.16) are then 
expressed as: 
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a r 


(3.2D 


r_ 
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6 6 • 6 
a a r 


r/6 


a 


6 

r 


6 r • 6 
r a 


<5 G 
r 6 6*6 
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•6 

a 


( 3 . 22 ) 


which are in the form of the single-input/single-output relation of 
Eq. (3.8) with cross-spectrum replaced by conditioned cross-spectrum, 
input-autospectrum replaced by conditioned input-autospectrum, and 
frequency response replaced by conditioned frequency-response. 

A generalization of the ordinary coherence function concept pre- 
sented above allows the quality of the conditioned frequency-response 
estimates to be assessed. The partial coherence function is defined as 
that fraction of the output that can be linearly related to a single 
input with the linear effects of all remaining inputs removed. It is 
interpreted in exactly the same way as the ordinary coherence functions 
for single-input/single-output systems described earlier Eq. (3.9), and 
has the same form. For the coupled system of Fig. 3.6, the partial 

coherence of yaw response to aileron with the linear effects of the 

2 

rudder inputs removed (y fi ) is calculated by 


|G 


6 r • 6 1 
a r 


'6 r * 6 
a r 


G 


6 6 *6 1 
a a r 


rr • 6 


(3.23) 


where G. . and G. . , are defined in Eqs. (3.17) and (3-18), 

& r*6 6 6 *6 

a r a a r 

respectively; and where 


G 

rr • 


6 


r 


e G (1 - y, ) 
rr 6 r 

r 


(3.24) 
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is termed the "conditioned output autospectrum" because it is the ordi- 
nary output autospectrum (G pr ) with the linear effect of partially 
correlated rudder excitation removed. The corresponding partial coher- 
ence of yaw-rate response to rudder with the linear effects of aileron 
2 

inputs removed (y. . ) can be obtained by interchanging 6„ and in 

6 r r*« a r a 

Eq. (3.23). 


Another useful metric is the multiple coherence function . This is 
the fraction of the output that may be linearly correlated to all of the 
inputs being considered. For the two- input case, the multiple coherence 

(y 6 6 :r J is given b y f 

a r * 


(r/5 a )0 « r 
a 


< r/ V G S r 

r 


'6 6 :r 
a r 


rr 


(3.25) 


The multiple coherence will be less than unity when input and output 
noise or nonlinearities are present, or when additional secondary con- 
trol inputs contribute to the aircraft response but are not included in 
the multi-input analysis. As before, the identified transfer functions 
will incorrectly appear decoupled when these additional secondary inputs 
are not taken into account. 

The two- input/s ingle-output spectral analysis calculations of 
Eq§~; (3. 15)-(3.25) may be performed using a computer program (TISOSA) 
developed by the author. The inputs to the program are the single- 
input/single-output spectral-response results of FRESPID, and a command 
file which contains the execution directives. The program outputs are 
the conditioned spectral -response plots and a data file. 


t Asterisk (*) denotes complex-conjugate. 
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3.3 Transfer-Function Modeling 


With the flight-test frequency responses in hand, the next step in 
the identification method of Fig. 1.3 is the derivation of transfer- 
function models. This step involves three important considerations 
which are treated in the next sections: 

1 . Order and form of transfer-function model 

2. Fitting range and relative weighting of the magnitude and phase 
curves 

3. Multi-degree-of-freedom models 

3.3.1 Selection of Model Order and Form 

Analytical transfer-function forms are selected for each input- 
output pair. These are based on configuration and flight condition 
considerations and must have a sufficient number of free parameters to 
adequately fit the identified frequency response over the frequency 
range of concern. In the hover flight condition, vehicle (rigid body) 
dynamics are dominated by the hovering cubic and decoupled heave and yaw 
modes. In wing-borne flight, the conventional longitudinal and lateral 
quartic equations dominate. Thus, transfer-function models which are 
appropriate for the hover flight condition are not necessarily appli- 
cable to forward-flight conditions. Obviously, if a model of high 
enough order is selected, the parameters can be adjusted to accommodate 
each flight condition. However, such transfer-function models no longer 
retain the physical significance of the classical lower-order parameters 
and are often not unique (Ref. 6). Also, higher-order transfer- 
functions models tend to be strongly tuned to the specific inputs which 
are used in the identification procedures (e.g., a frequency-sweep) and 
can be very poor predictors of responses to other test inputs (e.g., 
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step and pulse inputs) and responses at nearby flight conditions. 
Therefore, higher-order models are not desirable. The appropriate 
transfer-function order is the minimum order which can satisfactorily 
fit the frequency responses, with the upper limit taken as the physical 
order of the system. 

3.3.2 Transfer-Function Fitting Range and Relative Weighting 

The maximum frequency range usable for transfer-function fitting 
corresponds to the region of acceptable coherence. Often, however, the 
frequency range of suitable coherence extends beyond the range of appli- 
cability of the selected transfer-function models. For example, in the 
XV-15 cruise identification (Chap. 5), a simple short-period model is 
assumed, since it adequately characterizes the rigid-body pitch dynamics 
in the broad frequency range of 0.3-10 rad/sec. At the lowest input 
frequencies (0. 1-0.3 rad/sec), the vehicle response is dominated by the 
lightly damped phugoid dynamics; these dynamics are not very important 
for piloted handling-qualities related to attitude control tasks and so 
can be ignored for the present study. By restricting the fitting range 
to 0.3-10 rad/sec, where the short-period model is appropriate, the 
physical significance of the parameters is maintained. In other cases, 
structural resonances may be present at the high-frequency end of the 
range of good coherence. If a rigid body transfer-function model is 
assumed, the fitting range must be restricted to exclude the effects of 
such resonances (Ref. 6). 

The interactive computer program NAVFIT, originally developed at 
McDonnell Aircraft Company by Hodgkinson and Buckley (Ref. 37) and 
enhanced by the author to access the output file from the FRESPID pro- 
gram, fits the frequency-response data with a selected transfer-function 
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model. The numerator and denominator orders are arbitrary, and a time 
delay can be included to approximate the effects of unmodeled high- 
frequency dynamics. Least-squares fitting is completed using up to 
50 points spaced linearly across the logarithmic frequency scale in the 
desired frequency range. The relative weighting of the magnitude and 
phase errors is arbitrary, but a commonly selected value for handling- 
qualities applications is 1 dB magnitude error to 7 deg phase error^ 
(Ref. 17). The outputs from the NAVFIT program are closed-form 
transfer-function models and plotted frequency-response comparisons of 
the flight data and model. 

3-3.3 Fitting Multi-Degree of Freedom Frequency Responses 

A good review of the practical considerations involved in transfer- 
function fitting is given by Bischoff and Palmer (Ref. 15). They found 
that satisfactory identification of the transfer-function parameters 
generally requires the simultaneous consideration of at least two 
degrees-of-freedom. By constraining the transfer-function models to 
contain mutual denominator factors (as they do physically), the location 
of numerator factors becomes much better defined. In the cruise condi- 
tion, for example, the denominators of the longitudinal transfer func- 
tions relating elevator inputs to normal-acceleration and pitch-rate 
outputs (a z /6 e and q/6 e ) are the same. Simultaneously fitting these 
degrees-of-freedom maintains the commonality of the denominator factors 
(the short-period mode), and avoids the "galloping L a " problem (this 
refers to the problem of a poorly defined pitch-rate numerator location 


^Professor Stephen P. Boyd notes that this is the natural weighting that 
is obtained when the fit is performed on a real -imaginary plot of the 20 log.Q 
values of the complex transfer-function. 
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in single- input/s ingle-output q/6 g identification, Ref. 17). Two com- 
puter programs (LONFIT and LATFIT) were developed at McDonnell Douglas 
Aircraft Company by Givan, LaManna, and Hodgkinson (Refs. 38, 39) and 
modified by the author to allow multiple-degree-of-freedom fitting of 
the frequency-response output files from FRESPID. A new option incorpo- 
rated by the author into LONFIT allows the identification of the instan- 
taneous center-of-pitch rotation, x na (Sec. 5.2.2). 

eg 

3.4 Time-Domain Verification of Transfer-Function Models 

Time-domain verification involves driving the extracted transfer- 
function models with the flight data for inputs not used in the identi- 
fication procedure and comparing predicted and actual time histories. 
This analysis is necessary to check a number of important assumptions 
and steps in the frequency-domain method: 

Linearity assumption: The adequacy of linearized lower-order 

models in describing the time-domain dynamics of nonlinear systems 
is always a key question — especially for helicopters in hovering 
flight. Since the reference condition is zero velocity, even small 
translational-velocity perturbations can produce significant changes in 
the dynamic response characteristics. Lateral and longitudinal fre- 
quency sweeps can involve translational velocities of up to 40 ft/sec at 
low frequencies, so a unique lower-order linear model may not be achiev- 
able for the entire frequency range (Ref. 22). As discussed earlier, 
the extracted transfer-function model is actually a describing function 
that is strictly valid only for the input amplitudes which are used in 
the flight test experiment. However, since the motions encountered 
during frequency-sweep testing are typical of those in near-hover 
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maneuvering flight, the extracted describing functions are generally 
quite suitable. Time-domain verification checks this assumption. 

Transfer-function fitting mismatch: Sometimes, apparently small 

mismatches in the frequency-domain fit can produce surprisingly large 
discrepancies in time-domain correlations, especially when these mis- 
matches occur in the frequency range near the dominant vehicle modes. 

If the identified model poorly predicts the time-domain verification 
response, the transfer-function model order and structure is adjusted to 
improve the correlation. 

Overtuned models: Identification procedures tend to be tuned to 

the type of test inputs used (e.g., a frequency -sweep ) , and their accu- 
racy in predicting responses to other classes of inputs (e.g., steps 
and doublets) is often uncertain. The possibility of overly-tuned 
models increases as the order of the transfer function is increased 
(Ref. 40). By testing the transfer-function model with other input 
forms, overtuning is avoided. 

Identification of the open-loop frequency response from closed-loop 
flight data: Although a complete 90-sec frequency-sweep cannot be 

generated on the open-loop unstable XV- 15 aircraft in hover, shorter 
duration open-loop step and pulse inputs are possible. Comparisons of 
these SCAS-off responses with the extracted transfer-function models 
verifies that the true (unbiased) open-loop vehicle response, instead of 
the inverse feedback transfer-function (at worst), has been identified. 

3.4.1 Time-Domain Verification Method 

Time-domain verification may be done by converting the transfer 
functions to canonical form: 
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( 3 . 26 ) 


x(t) = Fx(t) + GUp(t)' 


y(t) = Cx(t) 


where 



and a Q , a 1f . . a n _ 1 are the denominator (ascending coefficients) 
and b Q , b^, . . . , b m are the numerator (ascending) coefficients of 
the transfer-function model being evaluated. 

The control input to the system u p is the low-pass filtered time- 
history of the surface deflection signal (e.g., 6 a in Fig. 3.7). The 
filter cutoff frequency is selected to eliminate spectral components of 
the test inputs (steps and pulses) which are above the frequency range 
of applicability of the derived transfer-function models (e.g., the 
upper bound of the least-squares fitting range, Sec. 3.3.2) (Ref. 6). 

State derivative biases and output reference values are included as 
unknown constants in the perturbation model of Eq. (3.19)> representing 
drifts and zero shifts, to allow a comparison with the flight responses: 
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(3.27) 


x(t) = Fx(t) + Gu„(t) + bx 

y(t) = cx(t) + y ref 

These bias and reference values are determined from a least-squares fit 
of the model response y(t) and the low-pass filtered flight data y p (t) 
[with x(0) = 0 to avoid overparameterizing the model]. The more 
common practice of simply subtracting the initial values of control and 
output from the time-histories can be used when the vehicle dynamics are 
stable, since drifts and reference errors decay with the time-constant 
of the transfer-function poles. However, when the identified transfer 
function is highly unstable, as it is in the pitch and roll dynamics of 
the XV-15 in hover, small errors in the biases and reference values 
quickly propagate and the responses diverge even for a highly accurate 
model. A computer program (TFTHISTORY) was developed by the author to 
estimate state-derivative biases and output reference values using a 
least-squares procedure (a simpler but still satisfactory alternative to 
a Kalman filter or smoother estimate), generate the model time-history, 
and compare the result with the filtered flight time-history data. 

This completes the discussion of the frequency -domain identifica- 
tion method. Identification results for the hover and cruise flight 
conditions are presented in Chapters 4 and 5. 
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Chapter 4 
HOVER DYNAMICS 


This Chapter discusses the identification and analysis of the XV- 15 
open-loop dynamics in hovering flight. The transfer functions of pri- 
mary interest for handling-qualities are associated with "inner-loop" 
regulation of the on-axis responses: pitch-rate response to elevator 

(q/6 e ), vertical-acceleration response to power lever (a 2 /« c ), yaw-rate 
response to rudder (r/« r ), and roll-rate response to aileron (p/« a ). 

The important coupling transfer function is yaw-rate response to aileron 
( r/6 ). To illustrate the identification method described in Chapter 3, 
the lateral/directional on-axis and coupled dynamics are discussed in 
detail. Flight-extracted frequency-responses are compared with piloted- 
simulation responses. The simulator responses are from the model 
states, so the drive dynamics of the vertical motion simulator (VMS) 
motion-base system are not included. Comparisons of the flight 
responses and the current nonreal-time simulation model (Ref. 3) are 
also included to show the substantial modeling improvements that have 
been made in the 3 years since the VMS simulation tests were conducted. 
Closed-form transfer-function models of the flight responses are 
obtained using least-squares fits. Finally step-response time-history 
matching is used to verify the predictive capability of the identified 
transfer functions. 

The results presented in this Chapter are considerably improved 
over the original hover flight condition results presented in Ref. 18. 
Subsequent experience on the XV-15 and related identification efforts 
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(Refs. 6, 21-24) led to refinements in the analysis methods and insight 
into techniques for extracting higher-quality results. 

4.1 Lateral/Directional Dynamics 

The frequency-sweeps in roll were the most difficult to execute 
because of significant interaxis coupling between lateral stick inputs 
and yaw response. Since the yaw SCAS employs a washed-out yaw-rate 
feedback, low-frequency yaw-attitude excursions required continuous 
pilot regulation during the lateral stick frequency-sweeps. However, a 
two-input analysis of the roll-axis response showed very little effect 
of the pedal inputs on the single-input aileron transfer function. This 
suggests that rudder inputs produce largely decoupled yawing motions; 
this is a typical assumption for hovering aircraft (Ref. 41). There- 
fore, the following analysis of the roll dynamics considers only the 
single-input/single-output transfer function p/S a . 

4.1.1 On- Ax is Responses 

The lateral-stick inputs for two good frequency-sweeps were concat- 
enated to form the composite input time-history of Fig. 4.1. The 20-sec 

FIRST SECOND 
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low-frequency cycles result in large lateral-velocity variations 
(10-20 knots). At the higher frequencies, stick deflections are larger, 
roll rates are smaller, and velocity perturbations are less than 
5 knots. Both runs are roughly 90 sec in duration, but the distribution 
of frequency content in each is seen to be quite different. The first 
frequency-sweep contains more mid-frequency power, whereas the second 
contains more low-frequency power . The input autospectrum for the 
concatenated time-history is shown in Fig. 4.2. Maximum input power is 



Fig. 4.2. Lateral-Stick Input Autospectrum ^ G 6 LAT « LAT ^’ 

generated at a frequency of u> = 0.4 rad/sec, which is slightly higher 
than that corresponding to the maximum period of T max = 20 sec. The 
input power remains strong out to w = 6 rad/sec, and then falls off 
rapidly at the higher frequencies. The frequency content which exists 
below the minimum even- input frequency (1/20 Hz) and above the maximum 
even- input frequency (4 Hz) (not shown in Fig. 4.2) is due to the 
various nonsinusoidal components which are apparent in the pilot input 
(Fig. 4.1). 

The corresponding aileron surface deflection time-history is shown 
in Fig. 4.3. This signal reflects inputs from the pilot (and shaped by 
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Fig. 4.3. Aileron Surface Deflections (6 a ) during the Lateral- 
Stick Frequency Sweeps, Reflecting Contributions From Pilot and SCAS 
Inputs. 

SCAS feed-forward elements) and from the SCAS feedbacks. Once again, 
the first frequency-sweep has more mid- and high-frequency power, and 
the second frequency-sweep has more low-frequency power. The input 
autospectrum for the concatenated surface deflection signal is shown in 
Fig. 4.4. The autospectrum peaks at to = 0.3 rad/sec (= low-frequency 



input) and is roughly constant up to 6 rad/sec. The fact that both the 
lateral-stick input autospectrum and surface-deflection input auto- 
spectrum are strong indicates that the surface deflection contains a 
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significant component from the pilot's stick inputs (thus avoiding the 
-1/Gj> problem discussed in Chap. 3). 

The roll-rate responses for the two lateral frequency-sweeps are 
shown in Fig. 4.5. The roll-rate amplitudes vary from 5-10 deg/sec, 
depending on input frequency. The output autospectrum (Fig. 4.6) peaks 



Fig. 4.5. Roll-Rate Response (p) during Lateral-Stick 
Frequency-Sweeps . 



Fig. 4.6. Roll-Rate Output Autospectrum (G p p). 


between 0.3 and 0.5 rad/sec, because of the low-frequency input cycles 
and the presence of the dominant roll modes in this frequency range. 
Although the input surface autospectrum remains roughly constant for 
w > 1 rad/sec, the output spectrum falls off rapidly because of the 
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rigid-body (inertial) response of the aircraft. The cross-spectrum 
between aileron and roll rate (Fig. 4.7) is quite similar in form to the 
roll-rate output autospectrum (Fig. 4.6) because the aileron input 
autospectrum (Fig. 4.4) is roughly constant with frequency. 



Fig. 4.7. Magnitude of Cross-Spectrum between Aileron Surface 


Inputs and Roll-Rate Response 



). 


The frequency response of p/6 a is shown in Fig. 4.8. Once again, 
the magnitude response peaks between 0.5 and 1 rad/sec, a result of the 
dominant roll modes which are in this frequency range. For frequencies 
greater than 1 rad/sec, the Bode plot shows the classic -20 dB/decade 
roll-off caused by the K/s (inertial) roll-rate response of the air- 
craft. The phase increase in the vicinity of the magnitude peak 
(to = 0.3-0. 6 rad/sec) indicates that the dominant roll dynamics are 
unstable. At higher frequencies, the phase is roughly -90 deg because 
of the inertial response and the small high-frequency lags. 

The coherence function for the roll-rate response identification is 
shown in Fig. 4.9. Excellent coherence is achieved in the frequency 
ranges of 0.2-0. 3 rad/sec and 0.6-9 rad/sec. The coherence function 
drops sharply in the frequency range in which the input and output 
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FLIGHT DATA 

SIMULATOR DATA 

NON REAL-TIME SIMULATION 



FREQUENCY, rad/sec 

Fig. 4.8. Roll-Rate Response to Aileron (p/-6 a ). (a) Magnitude; 

(b) Phase. 


FLIGHT DATA 

SIMULATOR DATA 



FREQUENCY, rad/sec 


Fig. 4.9. Coherence Function for Roll-Rate Response Identification 



autospectrum are both very strong and the transfer function peaks 
(0.3-0. 6 rad/sec). This indicates that nonlinear behavior is present in 
the low-frequency unstable dynamics. Current research by the author and 
Kaletka (Ref. 22) shows that a major source of this nonlinearity is 
rolling-moment dependency on higher powers of lateral velocity [L(v^), 
L(v^)]. For frequencies below 0.2 rad/sec and above 9.0 rad/sec, the 
coherence function drops rapidly because of the lack of excitation 
outside of the frequency range of dominant pilot inputs. 

The spectral results presented thus far are impressive in terms of 
their smoothness, resolution, and coherence. The chirp z-transform, 
described earlier, is a significant improvement over the standard 
(Cooley-Tukey) FFT methods, which were found to produce much poorer 
frequency-response resolution and associated erratic coherence-function 
behavior at low frequency. This is because resolution of the standard 
FFT is limited by the length of the window, whereas the chirp 
z-transform has arbitrary resolution (Sec. 3.2.2). In addition, the 
concatenation of good repeat runs is seen to produce a reasonably con- 
stant input (surface) autospectrum over the frequency range of interest, 
and provides a large number of sections for the averaging process. This 
ensures minimal random error and satisfactory accuracy of the transfer- 
function estimates, even when the magnitude of the coherence function is 
low (for example near u> = 0.4 rad/sec in Fig. 4.9). 

The hover frequency-sweep tests and analysis described above for 
the flight vehicle were repeated in the VMS. The extracted aileron-to- 
roll-rate open-loop transfer function is shown in Fig. 4.8 for compari- 
son with the flight results. The associated coherence function shown in 
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Fig. 4.9 reflects excellent simulator identification for the entire 
range of 0.1-10 rad/sec. 

The correlation of the VMS simulation and the flight data is very 
good for the frequency range 1.0-10.0 rad/sec, indicating that the roll- 
control-response sensitivity is accurately modeled. The low-frequency 
phase comparison suggests that the damping ratio of the unstable roll 
mode is slightly overestimated by the simulation model (more unstable), 
and the natural frequency is accurate. 

An interesting demonstration of the effects of nonlinearities in 
the low-frequency roll dynamics (p/6 a ) is presented in Fig. 4.10. This 
figure compares the simulator frequency-sweep result with the small- 
perturbation transfer function (obtained by numerically linearizing the 
VMS simulation model). The perturbation transfer-function matches the 
sweep results well in the high-frequency range (u» > 3.0 rad/sec) where 
motions are small; however, considerable deviation occurs in the low- 
frequency end, clearly a result of the large-perturbation nonlinear 
effects. (The large phase shift at low frequency is due to low- 
frequency, right-half plane poles and zeroes in the simulator perturba- 
tion transfer-function.) For forward-flight conditions, the differences 
were found to be much less significant since rotor nonlinearities are 
overshadowed by the linear aerodynamics of the wing and tail, and 
flight-path excursions about the steady-state trim condition are 
smaller. 

The fact that the simulation (and flight data) coherence function 
remains strong at low frequencies (Fig. 4.9) leads to the important 
conclusion that the dominant nonlinear effects can be well approximated 
by linear (but not small-perturbation) describing functions. By fitting 


63 


SIMULATOR DATA 

SIMULATOR PERTURBATION 
TRANSFER FUNCTION 




Fig. 4.10. Comparison of Small-Perturbation Transfer-Function and 
Frequency-Sweep Responses Obtained from the Simulation; Data Presented 
are for the Roll-Rate Response to Ailerons, p/-6 a - (a) Magnitude; 

(b) Phase. 

these describing functions with standard transfer-function forms, the 
physical significance of the small-perturbation model parameters is 
retained, and the numerical values of these parameters reflect the gross 
nonlinear effects. This "equivalent low-order modeling" approach has 
been useful in analyzing and predicting the handling qualities of many 
aircraft (Refs. 15, 17). 

The low-frequency errors in the simulation, which are shown 
in Fig. 4.8, were also reported during the pilot's qualitative 
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evaluation. Subsequent analysis of the real-time mathematical model 
indicated problems in the representation of rotor flapping for large 
lateral-velocity changes (such as occurs in the low-frequency range of 
the frequency-sweep). There is strong sensitivity of the numerically 
1 inear ized perturbation transfer-functions to the size of the perturba- 
tion increments. Corrections to the mathematical modeling in these two 
areas were included later in the nonreal-time XV-15 simulation. The 
comparison between the current nonreal-time mathematical model, the VMS 
original model, and the flight data is shown in Fig. 4.8; it indicates 
the significant improvements that have been achieved during the last 3 
years of development. 

The identified yaw-rate response to rudder inputs r/6 p is pre- 
sented in Fig. 4.11. The associated coherence function (Fig. 4.12) 
shows that the identification is accurate from 0. 1-3.0 rad/sec; above 
this range, the erratic coherence behavior suggests a less accurate 
estimate. 

The yaw-rate response follows the K/s characteristic in magnitude 
and phase for frequencies greater than 0.5 rad/sec, which means that 
rudder inputs generate constant yaw-rate accelerations. The fact that 
the phase does not rise above the -90 deg value except at very low 
frequencies suggests an associated low value of yaw damping (N r ) as 
compared to conventional helicopters. This is due to the lack of a tail 
rotor for the tilt-rotor configuration. 

A comparison of the simulator yaw-rate response with the flight- 
test data is presented in Figs. 4.11 and 4.12. As in the roll axis, the 
magnitude response generally matches very well for higher frequencies 
a) > 1.0 rad/sec. The parallel (vertical) shift in the magnitude curve 
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Fig. 4.11. Yaw-Rate Response to Rudder (r/6 r ). (a) Magnitude; 

(b) Phase. 
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FREQUENCY, rad/sec 


Fig. 4.12. Coherence Function for Yaw-Rate Response Identification 


indicates that the control sensitivity is overestimated by about 30 %. 
Subsequent analyses as a result of this finding indicate that ballast, 
which was added in the tail (for c.g. control), was not accounted for in 
the simulation model. This added inertia reduces the control sensitiv- 
ity and, therefore, could contribute to the observed high-frequency 
discrepancy. 

The yaw-mode damping is significantly underestimated in the simula- 
tor, as can be seen from the low-frequency phase drop (which actually 
indicates an unstable response!) compared with the (stable) flight 
response (Fig. 4.11b). This discrepancy was verified in step-response 
comparisons between the aircraft and the simulator. A likely cause of 
error in the yaw damping is the assumption of uniform lateral distribu- 
tion of the rotor inflow. Wing-on-rotor interference on the inboard 
half of the rotor disks probably causes a significant distortion in the 
shed wake and consequently in the lateral inflow distribution in hover- 
ing flight. The assumption of uniform inflow causes an underestimation 
in the longitudinal flapping response to longitudinal velocity changes 
(Ref. 42), which is the primary source of yaw damping for the hover 
flight condition. Even minor coupling to the (unstable) roll dynamics 
could then result in the observed unstable yaw simulation response. 

With the flight and simulation frequency-responses now identified 
and compared, the next step in the method of Chapter 3 is to derive 
lower-order equivalent transfer functions. Examination of the roll and 
yaw frequency-responses of Figs. 4.8 and 4.11 shows that they exhibit 
classic characteristics of hovering vehicles. The yaw response is 
decoupled and first-order in nature; the roll response is dominated at 
low-frequency by the unstable complex mode and has one excess pole at 
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high frequency. Thus, the yaw-rate response to rudder inputs is modeled 
by 


r(s) 

* r (s) 


N « 6 

r 


■t . s 


(s + 1/T y ) 


(4.1) 


where r(s)/6 r (s) is the Laplace-transformed yaw-rate response to 

rudder, in degrees per second per degree, and s is the Laplace- 

transform variable. The transfer-function high-frequency gain (N. ) is 

°r 

the yaw-control response sensitivity, that is, the initial yaw- 
acceleration response (deg/sec 2 ) to a 1-deg step input in rudder. The 
inverse yaw-mode time-constant (1/T y ) is the negative value of damping 
stability derivative (N r ). The effective time-delay (t^) accounts for 
phase distortion caused by unmodeled high-frequency dynamics . The 
numerical values for these parameters are obtained from the NAVFIT 
program (Sec. 3.3.2) using a fitting range of 0. 1-3.0 rad/sec: 


p 

= 0.619 deg/sec /deg-rudder 

p 

= 5.45 deg/sec^/in.-ped 

(4.2) 

1/T y = 0.102 rad/sec 
= 0.0210 sec 


As shown in Fig. 4.13, this model fits the flight data well at mid- 
and high-frequencies, but exhibits noticeable discrepancies at low 
frequencies, where the simple first-order time-constant approximation is 


68 



FLIGHT DATA 



FREQUENCY, rad/sec 

Fig. 4.13. Transfer-Function Model Identification for Yaw-Rate 
Response to Rudder (r/« r ). (a) Magnitude; (b) Phase. 

inadequate. The low-frequency mismatch is not very important, as is 
demonstrated in the model verification study (Sec. 4.3.1). The small 
effective time-delay (x^ = 21 msec) indicates that the unmodeled high- 
frequency dynamics (e.g., rotor and actuators) are not very important in 
the identification range (0. 1-3.0 rad/sec). 


69 


The standard roll-response transfer-function model for hovering 


vehicles is given in Ref. 41 as^ 


p(s) 

6 a (3) 


-T.S 

L. s( 1/T ) ( 1/T )e ♦ 

6 a *1 *2 

(1/T y )(1/T r )[c r ,«o r ] 


(4.3) 


where p(s)/6„(s) is the Laplace-transformed roll-rate response to 

a 

aileron surface deflection, deg/sec/deg-aileron ; L, is the aileron 

°a 

roll sensitivity; 1/T^ , 1/T^ are first-order, numerator inverse time- 
constants; 1/T , is the first-order roll mode inverse time-constant; c r » 
o> r are the second-order roll-mode damping and natural frequency, 
respectively; and is the effective time-delay. The decoupled yaw 

mode 1/T is fixed at its previously identified [Eq. (4.2)] value 
(1/Ty = 0.102 rad/sec). The remaining numerator and denominator parame- 
ters are varied to obtain a good fit in the range of 0. 2-9.0 rad/sec. 

The resulting roll transfer-function parameters are 


t Shorthand notation: [e, w] implies s 2 + 2cws + u> 2 , e = damping ratio; 

to = undamped natural frequency (rad/sec); and (1/T) implies s + (1/T), 
rad/sec. 
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= -3.71 deg/sec 2 /deg-aileron 

= 15.1 deg/sec 2 /in.-6 Lat 

1/T. = -0.107 rad/sec 

♦l 

1/T. =0.412 rad/sec 

♦2 

1/Ty = 0.102 rad/sec 

(4.4) 

1/T r = 1.23 rad/sec 

C r = -0.418 

a> r = 0.447 rad/sec 

t. = 0.0313 sec 
4 » 

As shown in Fig. 4.14, the transfer-function model is generally a very 
good representation of the aircraft roll response. However, a close 
examination of the model-fit and flight data in the frequency-range near 
the dominant mode (u> = 0.5 rad/sec) shows the inability of the linear 
model to simultaneously match both the magnitude and phase responses. 
Based on the magnitude response alone , a somewhat more unstable (higher 
negative damping ratio) mode is indicated; but, based on the phase 
response alone , a slightly less unstable (smaller negative damping 
ratio) model is indicated. This inconsistency reflects the significant 
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Fig. 4. 1H. Transfer-Function Model Identification for Roll-Rate 
Response to Ailerons (p/-6 a ). (a) Magnitude; (b) Phase. 

nonlinearities which are in this frequency range. The linear model of 
Eq. (4.4) is thus a compromise, and reflects the current relative 
weighting: 1 dB magnitude error: 7 deg phase error. Increasing the 

weighting of the phase errors resulted in a less unstable damping ratio 
(t ), better phase matching, and poorer magnitude matching as 
expected. The current weighting was retained since it is consistent 
with practice in the handling-qualities community (Ref. 37), and as 
shown in Sec. 4.3 it leads to transfer-function models with good predic- 
tive capability. The capability to specify relative magnitude and phase 
error weighting is an advantage of the frequency-domain approach. 
Time-domain identification methods inherently weight the phase errors 
much greater than the magnitude errors, since small shifts in phase 


72 


between the flight-data and the model-response generate large changes in 
the time-domain cost function (Ref. 22). 

Some important conclusions can be drawn from the transfer-function 
model parameters of Eq. (4.4), and the location of the poles and zeros 
in the complex plane (Fig. 4.15). First, the open-loop roll response 



Fig. 4.15. Locations of Roll-Response Transfer-Function Parameters 
in the s-Plane. 

dynamics are dominated by a very unstable oscillation in hover, with an 
associated time-to-double amplitude of 3.7 sec. Reference 43 specifies 
an unstable damping ratio limit of t = -0.22, for Level I handling 
qualities, which is shown on Fig. 4.15. The fact that the XV-15 roll 
mode [-0.418, 0.447] lies beyond this boundary suggests that the SCAS- 
off vehicle has Level II handling qualities (at best). This is cor- 
roborated by the SCAS-off handling-qualities evaluation discussed in 

p 

Sec. 2.1.3. The maximum roll acceleration capability of 1.27 rad/sec 
(based on the lateral-stick sensitivity and maximum deflection) provides 
a sufficient margin for SCAS-off roll-attitude control (AGARD specifica- 
tion, Ref. 43). 
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The second important conclusion is that the pole-zero combination 


(1/T. )/(1/T ) is at very low frequency and nearly cancels out 
<p 1 y 

(Fig. 4.15). This reveals that yaw coupling does not noticeably affect 
the roll-response characteristics, again a common assumption for hover- 
ing aircraft (Ref. 41). Third, the low-frequency numerator factor 

(1/T. ) associated with coupled lateral translation/yaw damping is also 
*1 

the limiting location of the dominant lateral translation mode when the 
inner attitude is tightly pilot-regulated (Ref. 44). The fact that this 
parameter is marginally unstable (time-to-double amplitude = 6.8 sec) 
suggests that lateral station-keeping will be difficult. This is typi- 
cal of unaugmented helicopters (Ref. 45) and is corroborated by the 
handling-qualities evaluation (Sec. 2.1.3). Fourth, the effective 
time-delay (t^) is small, again suggesting that unmodeled high-frequency 
dynamics are not significant in the frequency range of identification 


(0. 2-9.0 rad/sec) . 

4.1.2 Lateral/Directional Coupling 

As discussed previously, the dominant source of interaxis coupling 
for the hovering flight condition is yaw-rate response to aileron 
inputs . Identification of this coupling characteristic is important for 
assessing the pilot's work load in open-loop lateral station-keeping 
tasks . 

Clearly the pilot and SCAS feedbacks to the pedals will have a 
first-order effect on the identification of the open-loop yaw-rate 
response to aileron. Therefore, the application of the two- input/ 
single-output analysis method is essential. This case will be presented 
in detail, since it is an excellent example of how coupled bare-airframe 
characteristics can be extracted from a multi-loop system. 
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The yaw-rate response to the lateral-stick sweeps of Fig. 4.1 is 
shown in Fig. 4.16. Excursions generally do not exceed ±5 deg/sec. 
This is due to the combined effect of the pilot and stability- 
augmentation stabilization of the yaw degree-of-freedom during the 
lateral-stick sweep. In essence, the pilot is attempting to eliminate 
the yaw-rate response to aileron — that which is now being identified. 



Fig. 4.16. Yaw-Rate Response (r) during Lateral -Stick Frequency Sweeps. 


The autospectrum for aileron surface inputs (Fig. 4.17) is very 
flat across the entire frequency spectrum, as desired. The aileron 


G *a«a-«r 



Fig. 4.17. Aileron Surface Input Autospectrum. 
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autospectrum G (Fig. 4.17) is not the same as that in Fig. 4.4 
6 6 
3 . 3 

(although the time-history data are identical) because the optimum 
window length (Sec. 3.2.1) was found to be longer for the single-axis 
identification of p/« a than F° r the coupled-system identification. In 
this case, the output is the yaw-rate signal, which has an autospectrum 
that is not very flat (Fig. 4.18). The cross-spectrum between aileron 



Fig. 4.18. Yaw-Rate Output Autospectrum. 

input and the yaw-rate response (Fig. 4.19) is relatively constant in 
the high-frequency range, and peaks at low frequency (u> = 0.3 rad/sec), 
as does the output autospectrum (Fig. 4.18). 

The frequency response of aileron to yaw-rate, with rudder inputs 
present (r/« a :6 r ) is shown in Fig. 4.20. The frequency response is 
relatively flat in the mid-frequency range, indicating a fixed relation 
between yaw-rate and aileron when the (pilot and SCAS) yaw-axis stabili- 
zation is active. The ordinary coherence between yaw-rate and aileron, 
shown in Fig. 4.21, has a very low value over most of the frequency 
range. This indicates the presence of other inputs which significantly 
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Fig. 4.19. Magnitude of Cross-Spectrum between Aileron Surface 
Inputs and Yaw-Rate Response. 

r /V« r 

«■/* a ' 5 r 




Fig. 4.20. Yaw-Rate Response to Ailerons, (a) Magnitude; 


(b) Phase. 



Fig. 4.21. Coherence Function for Yaw-Rate Response to Aileron 
2 

Inputs ( Yg r ). 

a 

contribute to the observed yaw-rate output signal; clearly, these are 
the rudder inputs, which are being supplied by the pilot and the yaw- 
stabilization system. 

Rudder surface inputs during the lateral frequency-sweep are shown 
in Fig. 4.22. The deflections are generally within ±5 deg, with most of 
the low-frequency activity coming from pilot inputs. The coherence 



Fig. 4.22. Rudder Surface Deflections (« r ) during the Lateral- 
Stick Frequency Sweeps. 
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function for the cross-correlation between rudder and aileron surface 
inputs presented in Fig. 4.23 is very low over most of the frequency 



Fig. 4.23. Coherence Function for Cross-Correlation between 

2 

Aileron and Rudder Surface Inputs (y. ). 

o 6 
a r 

range. This shows that the inputs are not fully correlated, the key 
requirement for two-input/single-output identification (Sec. 3.2.5). 

For other flight conditions in which similar analyses were conducted, 
the cross-coherence between rudder and aileron inputs was generally much 
greater because pilot-stabilization requirements in the secondary axis 
were less severe. This increased cross-coherence caused an overall 
degradation in the cross-coupling identification for the forward-flight 
cases. Therefore, when possible, the frequency-sweeps should be con- 
ducted without engaging the stability-augmentation systems in the sec- 
ondary degrees-of-freedom. 

The conditioned aileron input autospectrum (G ), that is, 

06 *6 
a a r 

with the effect of rudder inputs eliminated, is presented in 
Fig. 4.17. The comparison of results in this figure shows that the 
input autospectrum is not substantially changed. This is because the 
cross-coherence between inputs is small. A similar characteristic is 
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apparent for the conditioned output autospectrum ) shown in 

r 

Fig. 4.18. The conditioned cross-spectrum between aileron and yaw-rate 

(G ) presented in Fig. 4.19 is generally more flat than before, a 

Sr- 6 
a r 

result of the elimination of the yaw-rate feedback effect. 

The conditioned frequency-response of aileron to yaw-rate (r/6 *6 ) 

a r 

is shown in Fig. 4.20. Notice that the shape of the magnitude and phase 

responses is significantly different than it was before the rudder 

inputs were taken into account. The response from 1.0-10.0 rad/sec now 

exhibits the familiar -20 dB/decade roll-off. This suggests that the 

interaxis coupling is dominated by the aileron-to-yawing moment 

control-coupling derivative (N, ), as expected. The low-frequency 

°a 

dynamics are unstable, as seen earlier in the roll response. 

2 

The partial coherence function y . . of Fig. 4.24 shows that 

6 r-6 
a r 

satisfactory identification of the coupling dynamics has been achieved 



Fig. 4.24. Partial Coherence Function for Yaw-Rate Response to 
2 

Aileron Inputs {y . . ) . 

O P * O 

a r 

in the frequency ranges of 0. 1-0.4 rad/sec and 0.75-7.0 rad/sec. The 

p 

multiple coherence plot y^ of Fig. 4.25 indicates that the com- 

a°r ' 

bined linear effects of aileron and rudder account for most of the 
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Fig. 4.25. Multiple Coherence Function for Yaw-Rate Response to 

p 

Combined Aileron and Rudder Surface Inputs (y ). 

^ 6 6 :r 

a r 

yaw-rate response at high frequencies with a drop-off at lower frequen- 
cies because of nonlinearities and turbulence, as was seen previously in 
the single-input, roll-response analysis. 

The following transfer-function model adequately characterizes the 
aileron-to-yaw rate coupling: 


~ T , s 
*2 


r(s) 

6 a (s) 


V^/yt v% ]e 

(1/T y )(1/T r )[c r ,« r ] 


(4.5) 


where r(s)/6 a (s) is the Laplace-transformed yaw-rate response to 

aileron surface deflection, deg/sec/deg-aileron; N, is the aileron yaw 

°a 

sensitivity; 1/T^ is the first-order numerator inverse time-constant; 

and o)^ are the second-order numerator damping and natural frequency, 
respectively; and is the effective time-delay. The denominator 

factors are the previously identified lateral modes [Eq. (4.4)]. These 
denominator factors are held constant, and the remaining parameters are 


81 


iteratively varied using NAVFIT to fit the data in the frequency range 
0. 1-7.0 rad/sec. The resulting transfer-function parameters are 


N s = 0.344 deg/sec^/deg-aileron 
a 

= -1.40 deg/s^/in.-6^ a j. 

1/T^ = -0.345 rad/sec 

= 0.868 

= 0.487 rad/sec (4.6) 

1/T = 0.102 rad/sec 

1/T p = 1.23 rad/sec 
t r = -0.418 
u> r = 0.447 rad/sec 


T 


4>2 


= 0.00900 sec 


The yaw-rate responses to both aileron and rudder are dominated by 
the "K/s characteristic" in the pilot's operating (crossover) frequency 
region for attitude control [generally taken as 1. 0-3.0 rad/sec 
(Ref. 41)]. Therefore, the pilot cross-feed (Gq F ) that is needed to 
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maintain a fixed heading during lateral manuevering is approximately a 
constant: 

. N * 

G CF = = deg-rud/deg-ail = 0.489 %-ped/^-lat (4.7) 

6 

r 

This shows that pedal inputs must be about 50% of the lateral stick 
inputs to maintain a constant heading in the piloted crossover frequency 
region. This level of coupling is quite severe, and probably contrib- 
utes to the poor SCAS-off handling-qualities in lateral translation 
tasks (Sec. 2.1.3). 

The comparison between the transfer-function model fit and the 
flight-test results for this case is shown in Fig. 4.26. The magnitude 
response correlates well with somewhat poorer agreement for the phase 
response. However, the nature of the coupling characteristics is well 
represented. 

In the simulated lateral-stick sweeps, pilot regulation of yaw 
excursions was not very aggressive. Therefore, the rudder inputs were 
dominated by the SCAS, causing a strong correlation between the rudder 
and aileron responses. This precluded the application of the two-input 
method. Instead, the flight-test results are compared in Fig. 4.26 with 
the real-time mathematical model (numerical) perturbation transfer- 
function. As expected, the high-frequency correlation is excellent, 
since it depends only on the cross-coupling moment derivative . The 

cl 

agreement at low frequency is much less satisfactory because nonlinear 
yaw and roll dynamics are not well represented by the small perturbation 
transfer-function ( Sec . 4.1.1). 


83 


FLIGHT DATA 

SIMULATOR PERTURBATION 

TRANSFER-FUNCTION 




Fig. 4.26. Comparison of Flight-Test and Simulator Results for 
Yaw-Rate Response to Aileron (r/$ a ). (a) Magnitude; (b) Phase. 


4.2 Longitudinal Dynamics 

Power lever inputs produce essentially decoupled first-order 
vertical-acceleration (heave) responses. So, these can be considered 
independently from the pitch responses. The vertical acceleration 
response to power-lever inputs (a z /6 c ) is shown in Fig. 4.27. For 
low-frequency inputs, the first-order heave mode dominates the response, 
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Fig. 4.27. Vertical-Acceleration Response to Power Lever 
(-a z /6 c ). (a) Magnitude; (b) Phase. 


causing the phase rise seen in Fig. 4.27b. The relatively slight change 
in the magnitude response (Fig. 4.27a) suggests that this mode is at a 
very low frequency, as is verified later. In the frequency range of 
1.0-10.0 rad/sec, the phase curve is nearly constant, so the aircraft 
vertical-acceleration response to power-lever inputs is instantaneous, 
and the high-frequency engine-governor dynamics are not significant. 

The coherence function is shown in Fig. 4.28, and indicates an excellent 
identification of the vertical response characteristics in the frequency 
range of 0. 1-3.0 rad/sec. 

The simulation-response characteristics are also shown in 
Figs. 4.27 and 4.28. In general, the agreement is good in the form 
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Fig. 4.28. Coherence Function for Vertical-Acceleration Response 

p 

Identification (y„ ). 

6 a 
c z 

of the amplitude and phase characteristics for frequencies 
o) > 0.3 rad/sec, suggesting that the modeling of initial thrust response 
to collective inputs, and governor dynamics, is satisfactory. The 
constant amplitude error (vertical shift) may be due to a mismatch in 
the assumed vehicle mass or rotor lift curve slope. The low-frequency 
agreement is not as good, which appears to be due to an underestimate of 
the heave-mode time constant by the simulation. 

The vertical-velocity response (w = a z /s) to collective is identi- 
cal in form to the yaw-rate response to rudder presented earlier. In 
fact, these two decoupled degrees-of-freedom are directly analogous. 
Typically, they both have very long time-constants for tilt-rotor-class 
VTOL vehicles, in contrast to single-rotor helicopters, so that their 
response is basically an acceleration; the control sensitivity is a 
function of the applied control force (or moment) normalized by the 
respective mass characteristic. A first-order model fit to the 
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vertical-acceleration response in the frequency range 0. 1-3.0 rad/sec 
yields 


a (s) 


-0.00980s e 


-0.00740s 


(s + 0.105) 


g/% 


(4.8) 


where the heave-mode time-constant has about the same value as the 
yaw-mode time-constant presented earlier, and the time delay is again 
small. The control sensitivity (-0.0098 g /% ) is equivalent to 
-0.1 g/in. of power lever; this is adequate, based on the data of 
Ref. 43. The very low value of the flight-extracted inverse time- 
constant (1/T h = 0.105 rad/sec) relative to typical hovering helicopters 
is due to the high disc-loading of the tilt-rotor configuration. 

Although the frequency-sweeps in roll were more difficult to per- 
form than those in pitch, the identification of the pitch results proved 
to be the most difficult. Stick inputs tended to be very small at low 
frequency, which resulted in poor uncorrelated input power. Of the 
three longitudinal frequency-sweeps which were conducted (and used in 
the original analysis of Ref. 18), only one run subsequently proved to 
be of satisfactory quality (based on input autospectrum and coherence) 
for the pitch-response identification in the mid- and low-frequency 
range. In the high frequency range (u> > 2 rad/sec), all sweeps are of 
good spectral quality, so concatenation of the time-histories is useful 
to obtain minimum variance results there. However, to achieve good 
identification over the entire frequency range only the single 
frequency-sweep is used. The possibility of combining spectral results 
which are optimized for particular frequency-ranges has not been fully 


87 


investigated, but may enhance the overall quality of the identification 
results. 

The elevator to pitch-rate frequency response is shown in 
Fig. 4.29. As before in the roll axis, the response peak is due to the 
dominant pitch modes which are at low frequency. For frequencies 
oj > 2 rad/sec, the magnitude response rolls off with the K/s rigid- 
body characteristic, as expected. The phase rise near the low-frequency 
response peak reflects an unstable pitch-mode associated with the 
speed/attitude divergence, as in the roll case. At high-frequency, the 
pitch-rate and elevator traces are just over 90 deg out of phase because 

of the rigid-body response and the small high-frequency lags. 

FLIGHT DATA 

SIMULATOR DATA 

MODEL 




Fig. 4.29. Pitch-Rate Response to Elevator (q/-6 e ). 
(a) Magnitude; (b) Phase. 
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The coherence function shown in Fig. 4.30 indicates a strong iden- 
tification in the frequency range 0. 2-6.0 rad/sec. The drop in coher- 
ence between 0. 5-1.0 rad/sec suggests that nonlinearities are important 
in low-frequency dynamics, as would be expected from the similar roll- 
response behavior. 

FLIGHT DATA 

SIMULATOR DATA 



Fig. 4.30. Coherence Function for Pitch-Rate Response 
2 

Identification (y r )• 

6 e q 

The simulation pitch-response characteristics are also shown in 

Figs. 4.29 and 4.30. The parallel (vertical) shift in the magnitude 

plots at high frequency suggests differences in the high-frequency gain 

parameter, M, . As in the yaw-case, the discrepancy (about 25 %) is 
°e 

probably a result of ballast in the tail of the aircraft, which was not 
included in the simulation model. Considerable differences between the 
simulation and flight-extracted responses are apparent in the frequency 
range of 0.3-0. 8 rad/sec. The simulation model appears to underestimate 
both the frequency and the unstable damping ratio of this mode. Com- 
parisons of time-domain identification results (Ref. 22) with the 
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nonreal-time simulation indicates that "speed-stability derivative" 

M u is underestimated in the simulator by about 70% 

t (Mu> FLIGHT = °-°°622 versus (^SIMULATION = °- 00200 rad/sec 2 /ft/sec ] 
which would contribute to the observed frequency and damping 
discrepancies (Ref. 42). A probable source of this error in the M u 
derivative is the rotor longitudinal flapping angle predictions based on 
the uniform lateral inflow approximation, as was also indicated in the 
yaw-rate response (Sec. 4.1.1). 

The transfer-function model for pitch-rate response to elevator has 
the same form as the roll-response model of Eq. (4.3): 


M « s(1/T )(1/T ) 6 

e e 1 e 2 


-v 


q(s) 

ys> - <i/T h ><i/T p )u p ,. p ] 


(4.9) 


where q(s)/S g (s) is the Laplace-transformed pitch-rate response to 
elevator surface deflection, deg/sec/deg-elevator ; M g is the elevator 
pitch sensitivity; 1/T 0 ^, 1/T 0 are first-order numerator inverse 
time-constants; 1/Tp is first-order pitch-mode inverse time-constant; 
Gp and u>p are the second-order pitch-mode damping and natural fre- 
quency, respectively; and t 0 is the effective time-delay. As before, 
the decoupled mode, (in this case 1/T h ), is fixed at its previously 
identified value (1/T^ = 0.105 rad/sec). The remaining numerator and 
denominator parameters are varied to obtain a good fit in the frequency 
range 0. 2-5.0 rad/sec: 
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= -2.66 deg/sec 2 /deg-elevator 
= -13.2 deg/sec 2 /in.-6 Lon 

1/T„ = -0.271 rad/sec 

e 1 

1/T„ = 0.508 rad/sec 

0 2 

1/T h = 0.105 rad/sec 

(4.10) 

1/Tp = 1.32 rad/sec 
C p = -0.463 
Up = 0.579 rad/sec 
t q = 0.0656 sec 

A comparison of the lower-order model with the flight-test results is 
shown in Fig. 4.29. Overall, the model is a good approximation to the 
flight data in the fitting range. However, as in the roll case, the 
presence of nonlinearities in the low-frequency unstable dynamics pre- 
cludes an entirely representative fit in both magnitude and phase near 
the response peak. 

Once again, some important conclusions can be drawn from the 
lower-order transfer-function model. First, the open-loop (bare- 
airframe) configuration is very unstable in the pitch axis, with nearly 
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the same damping ratio and natural frequency as in the roll axis. As 
before, poor SCAS-off handling-qualities are indicated and have been 
reported in piloted evaluations (Sec. 2.1.3). Second, the pole-zero 


pair associated with heave motions ( l/Tg^J/C 1/T h ) is widely separated, 

implying that significant heave-pitch coupling results from pitch- 

control inputs. This separated pole-zero pair precludes the use of the 

standard third-order, pitch transfer-function model for hover. Third, 

the numerator factor (1/T Q ) associated with the translational mode 

e 1 

dynamics for constrained attitude control (Ref. 44), is very unstable 
(time-to-double amplitude = 2.6 sec), and probably contributes to the 
poor translational-flight handling-qualities (Sec. 2.1.3). Fourth, the 
effective time-delay (t 0 ) is very small, as it is in the roll axis, 
suggesting that unmodeled high-frequency dynamics are not significant in 
the frequency range of identification (0. 2-5.0 rad/sec). 


4.3 Verification of Open-Loop Transfer-Function Models for Hover 

Verification of all identified open-loop transfer functions should 
be completed with open-loop flight-test data . However, the duration and 
size of the verification test-inputs is severely restricted by the 
inherent stability and response characteristics for each axis and flight 
condition. The decoupled stable first-order nature of the open-loop yaw 
and heave responses allows large and very long duration step-inputs 
(about 10 sec) to be executed. Then, the modes with frequencies of 
co > 0.1 rad/sec contribute at least one time-constant of their inherent 
response to the overall behavior, which checks the yaw and heave trans- 
fer-functions [Eq. (4.2) and (4.8), respectively] over the entire iden- 
tified frequency range. However, as previously discussed, the open-loop 
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pitch and roll dynamics are dominated by strong instabilities, with 
time-to-double amplitudes of about 3 sec. This severely restricts the 
allowable size and duration of the verification test inputs in these 
axes. 

To reduce the level of coupled divergence, verification of the 
open-loop pitch response is completed with the roll and yaw SC AS axes 
engaged . Similarly, verification of the open-loop roll response is 
completed with the pitch and yaw SCAS axes engaged . Even so, maximum 
achieveable step-input durations are 1-2 sec, which allows the dominant 
modes with frequencies of about w = 0.5 rad/sec to experience only 
one time constant. Therefore, the verification step- inputs do not give 
a complete test of model fidelity at very low frequencies. However, 
ensuring model accuracy at mid- and high-frequency is of principal 
importance for handling-qualities (and control) applications since 
the pilot is mostly concerned with the response in the first few 
seconds. SCAS-on step responses are often useful for verification when 
quiescent initial test conditions cannot be achieved in the open-loop 
configuration. 

4.3.1 Roll Response to Aileron Inputs 

The (filtered) roll-rate response of the aircraft to a step aileron 
input with the roll SCAS disengaged is shown in Fig. 4.31. The input is 
held for 2 sec and the peak response is 10 deg/sec, which is roughly the 
same as was achieved during the lateral frequency-sweeps (Fig. 4.5). 

When the transfer-function model of Eq. (4.4) is driven with the same 
(filtered) input flight data, the response co-plotted in the dashed 
curve of Fig. 4.31 is obtained. The model and flight data responses 
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Fig. 4.31. Roll-Response Model Verification Using SCAS-off Flight 
Data, (a) Aileron Input; (b) Roll-Rate. 

generally compare very well, although there is a slight underestimation 
of the peak roll rate, probably a result of nonlinear effects. 

4.3.2 Yaw Response to Rudder Inputs 

The yaw-rate responses of the aircraft and model Eq. (4.2) to an 
open-loop step-rudder input in hover are shown in Fig. 4.32. The con- 
trol deflection represents about 40% of the maximum control power, a 
fairly large amplitude input. The maximum yaw-rate response is about 
35 deg/sec, which is 40£ greater than the peak values obtained in the 
frequency sweeps. The model and flight responses generally compare very 
well over the entire run — even for this fairly large input amplitude 
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Fig. 4.32. Yaw-Response Model Verification Using SCAS-off Flight 
Data, (a) Rudder Input; (b) Yaw-Rate. 

case. The small discrepancies are probably associated with the mid- 
frequency mismatch between the flight-extracted frequency response and 
transfer-function fit (Fig. 4.13), or large amplitude nonlinear effects. 

4.3.3 Yaw Response to Aileron Inputs 

The verification of the identified model for yaw-rate response to 
aileron [Eq. (4.6)] is shown in Fig. 4.32. The SCAS configuration is 
roll, SCAS-off; and yaw, SCAS-on. Therefore, some rudder activity is 
also present (Fig. 4.33b). The roll-rate match is very good 
(Fig. 4.33c) as was seen in the previous comparison (Fig. 4.31). The 
observed yaw-rate (Fig. 4.33d) is composed of responses to rudder and 
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aileron inputs in about equal proportions. The total yaw-rate response 
has the proper form although the peak is underestimated; this is prob- 
ably a result of nonlinearities in the coupled response characteristics 
(or a gust), as well as the slight underestimation in the yaw-rate 
response to rudder (Fig. 4.32). 

4.3.4 Vertical Acceleration Response to Collective 

The verification data for the step inputs in the collective axes 
are very noisy. Even so, the correlation of the model and flight data 
(Fig. 4.34) indicates that a proper identification has been achieved. 
The comparison is better at the beginning of the run where the inputs 
are smooth as compared to the end of the run where the inputs are more 
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Fig. 4.34. Vertical-Acceleration Response Model Verification, 
(a) Power Lever Input; (b) Vertical Acceleration. 
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sharp. This is because the sharp change in collective (rotor thrust) 

excites the rotor inflow dynamics which are not included in the simple 

first-order model (Ref. 6). 

4.3.5 Pitch-Rate Response to Elevator 

The aircraft and model responses to a 2-sec step-input in elevator 
with the pitch SCAS disengaged are shown in Figs. 4.35. The relatively 

FLIGHT DATA 

MODEL 




Fig. 4.35. Pitch-Response Model Verification Using SCAS-off Flight 
Data, (a) Elevator Input; (b) Pitch-Rate. 

unsteady initial conditions cause the shift between the flight data and 
transfer-function models at the beginning of the run (t = 6-8 sec), but 
the comparison improves considerably during the large amplitude recovery 
phase (t = 8-15 sec). Somewhat better initial conditions are achieved 
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with the pitch-SCAS engaged. The elevator surface-deflection response 
to a step input of longitudinal stick is a doublet, as seen in 
Fig. 4.36a. As expected from the control system configuration 
(Chap. 2), a steady pitch-rate is achieved within 0.5 sec of the longi- 
tudinal stick input. The model prediction is seen to be very good for 
the input and recovery (t = 9-13 sec). 

FLIGHT DATA 




Fig. 4.36. Pitch-Response Model Verification Using SCAS-on Flight 
Data, (a) Elevator Input; (b) Pitch-Rate. 
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The verification results confirm the following: 

1. The hover dynamics are essentially linear. 

2. Mismatches between the extracted frequency-responses and 
transfer-function model fits are not significant, except perhaps in the 
yaw-rate response to rudder. 

3. The derived transfer-function models are not overtuned or 
over-parameterized. 

4. The unbiased open-loop vehicle response has been identified 
from the closed-loop flight data. 

The derived transfer-function models for hover are summarized in 
Table 4.1. The verification results show the utility of these fairly 
simple models in representing the transient responses to relatively 
large and varied inputs. Therefore, the models presented in Table 4.1 
should be used for future XV-15 control system and handling-qualities 
analyses for the hover flight condition. 

The next chapter discusses the dynamics for the cruise flight 
condition. 
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TABLE 4.1 Identified Transfer-Function Models 


for Hover 
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Chapter 5 
CRUISE DYNAMICS 


This chapter discusses the open-loop dynamics of the XV- 15 for the 
cruise flight condition (flight condition 4, Sec. 2.2). 

As in hovering flight, the longitudinal and lateral flight dynamics 
are essentially fully decoupled. The primary longitudinal bare-airframe 
transfer functions of interest for "inner-loop" regulation are pitch- 
rate and normal-acceleration (at the center-of-gravity ) responses to 
elevator deflection (q/6 g ) and (a z /6 e ) respectively. The important 
lateral-directional transfer functions are roll-rate response to aileron 
(p/6„) and sideslip (at the center-of-gravity) response to rudder 

3. 1 

IV'r 1 ’ 

Frequency-sweep tests were executed much more easily in the cruise 
flight condition than in the hover condition. Cruise dynamics are 
inherently very stable, and transient motion is small relative to the 
large steady airspeed (V^ = 170 knots). This allowed good flight data 
to be rapidly obtained without requiring a large number of repeti- 
tions. Also, high-quality spectral identification was achieved without 
considerable adjusting of the various processing options — as was neces- 
sary in hover. Only two frequency-sweeps were needed in each axis 
because the spectral quality of all of the individual runs was satisfac- 
tory, and none needed to be discarded (as was required in the hover 
case). The chirp z-transform algorithm yielded excellent frequency- 
response identification in both axes over a wide frequency range. 

Frequency-sweep data from the (VMS) piloted simulation were not 
available for the (V^ = 170 knots) cruise flight condition, so a direct 
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comparison was not possible. However, the nonreal-time mathematical 
model (Ref. 3) for this condition is essentially identical to the real- 
time piloted-simulation (Ref. 2). Since the relative amplitude of 
motion about the high equilibrium flight speed is generally small, and 
the dynamics in the fixed-wing configuration are basically linear, the 
perturbation transfer-functions obtained from the nonreal-time version 
are satisfactory for mathematical model validation in this flight 
condition. 

5.1 Lateral/Directional Dynamics 

5.1.1 Frequency-Response Identification 

The lateral stick inputs for two concatenated frequency-sweeps in 
the hover flight condition are shown in Fig. 5.1a. The input form is 
seen to be much more regular than in hover (compare with Fig. 4.1). The 
two runs are quite similar, with maximum amplitudes of about 20 % of the 
full control authority. Since these tests were conducted with the 
lateral/directional SCAS axes disengaged , the surface deflection (exci- 
tation input) is related to the cockpit signal by a constant (mechanical 
compliance and hydraulic actuators contribute less than 10 msec of 
effective time delay). The roll-rate response for these two frequency- 
sweeps is shown in Fig. 5.1b. The maximum response amplitude of 
10 deg/sec is similar to that experienced in hover (Fig. 4.5). 
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TIME, sec 

Fig. 5.1. Two Lateral Stick Frequency-Sweeps ( 5 lat ^ Cruise. 

(a) Lateral Stick Inputs; (b) Roll-Rate. 

The (open-loop) roll-rate response to ailerons, p/6 a , is shown in 
Fig. 5.2. This is a classical first-order response in which aileron 
inputs produce a constant roll acceleration at high frequency and a 
constant roll rate at low frequency. The corner frequency is at about 
to = 1 rad/sec, with an associated phase lag of roughly ♦ = -45 deg, as 
expected. The dominant time-constant is thus about 1 sec, implying 
about 2-3 sec to reach a steady-state roll rate. Higher-order dynamics 
are not significant in this response since the phase lag is nearly 
constant for frequencies above 3 rad/sec. The coherence function shown 
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Fig. 5.2. Roll-Rate Response to Aileron (p/-S a ). (a) Magnitude; 

(b) Phase. 

in Fig. 5.3 is very strong over the frequency range 0. 1-8.5 rad/sec, 
indicating excellent identification of the roll-rate dynamics. 

The perturbation transfer-function response from the nonreal-time 
simulation is shown with the flight data in Fig. 5.2. The magnitude 
response comparison is excellent at high and low frequencies, and shows 
only a slight error at mid-frequencies. The phase comparison 
(Fig. 5.2b) indicates that the simulation model overpredicts the damping 
ratio of the dominant roll-response mode (<j = 1 rad/sec). The phase 
error at high frequency is due to the omission of the high-frequency 
actuator dynamics in the simulation model. These various small discrep- 
ancies in Fig. 5.2 are insignificant because it is not possible to 
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Fig. 5.3. Coherence Function (y ) for Roll-Rate Response 

a^ 

Identification. 

perfectly match the test conditions of the aircraft with the simula- 
tion. Clearly, though, the simulation response characteristics are a 
very good representation of the actual vehicle response. 

Aerodynamic sideslip measurements are obtained from a sideslip 
indicator located about 18 ft ahead of the aircraft center-of-gravity 
(c.g.). The sideslip at the c.g. (B C g) is calculated by correcting the 
measured signal for position error based on yaw-rate and airspeed (no 
flow-distortion corrections were made). Also, no corrections are made 
for sensor dynamics, because these are not felt to be significant within 
the identification bandwidth. 

The sideslip response (at the c.g.) to rudder inputs, 8 C g/^ r > i- s 
shown in Fig. 5.4. The response is characterized by a lightly damped 
second-order mode with a frequency of about 1.6 rad/sec. High-frequency 
rudder inputs yield a constant sideslip acceleration with no appreciable 
time-delay, and low-frequency rudder inputs yield a constant sideslip 
angle. The coherence function (Fig. 5.5) is strong over the frequency 
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range 0. 1-3.5 rad/sec, falling off sharply at higher frequencies. This 
rapid and relatively low-frequency decline of the coherence function 
results because the sideslip response falls off at a rate of 
-40 dB/decade, in contrast to the -20 dB/decade roll-off for the other 
response variables. However, since the yaw mode of interest has a 
natural frequency of about 1.6 rad/sec, the falling coherence for fre- 
quencies greater than 3.5 rad/sec is not a severe limitation. 

The simulation frequency response of 6 C g/6 r is compared with the 
flight data in Fig. 5.4. As in the roll case, the comparison is gener- 
ally excellent. The parallel (vertical) shift in the magnitude curves, 
suggests that the rudder yaw-response sensitivity is overestimated in 
the mathematical model by about 30$. This matches the rudder-response 
sensitivity discrepancy in hover and so further indicates a yaw inertia 
error. The slightly more shallow phase roll-off in the simulation 
response at high frequency is again due to the omission of actuator 
characteristics and perhaps a small overestimation in the damping ratio 
of the Dutch roll mode. Once again, however, these differences are not 
considered to be significant and could be attributed to differences in 
flight and loading conditions between the aircraft and the simulation 
model . 

5.1.2 Transfer-Function Models 

Equivalent system fitting for fixed-wing aircraft using decoupled 
lateral/directional models has been considered in detail by Bischoff and 
Palmer (Ref. 15). With the decoupled model, the responses are fit 
independently, using a first-order roll rate and a second-order sideslip 
transfer function. The coupled model approach is based on simultaneous 
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fitting of the roll and sideslip responses to aileron and rudder, 
respectively, i.e., we obtain the following fourth-order transfer 
functions: 


L. r i V 

p(s) v [ v % le 

5 a < = > - (1/T a )(1/T r )U d , . d ] 


(5.1) 


where p(s)/6 a (s) is the Laplace-transformed roll-rate response to 

aileron surface deflection (deg/sec/deg-aileron) , L g is the aileron 

s 

roll sensitivity, and are the equivalent Dutch roll-mode damping 
and natural frequency, respectively, and is the effective time 

delay, and 
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where B C g(s)/S r (s) is the Laplace-transformed sideslip response to 
rudder surface deflection (deg/deg-rudder) ; Y g is the rudder sideslip 
sensitivity; 1/Tg^, 1/Tg^, and 1/Tg are the first-order numerator 
inverse time constants; and the denominator parameters are identical to 
those of Eq. (5.1). The effective time delay for the sideslip response 
is Tg. The simultaneous fitting approach is consistent with that used 
in the longitudinal axis and allows the identification of the 
"coupling effect" in the roll response [Eq. (5.1)], which is important 
in handling-qualities assessments. 

The coherence function results of Figs. 5.3 and 5.5 show a satis- 
factory identification of both roll and sideslip responses in the 
frequency range 0. 1-3.5 rad/sec. Based on this fitting range, the 
computer program LATFIT is used to obtain the parameters of Eqs. (5.1) 
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satisfactory over the wider range of 0. 1-8.5 rad/sec, this degree of 
freedom is refitted alone, holding the denominator factors constant at 
the values obtained from the simultaneous solution. This procedure 
optimizes the values of the high-frequency gain and numerator parameters 
in the roll transfer function. The final results for the parameters of 
Eqs. (5.1) and (5.2) are 


2_/ \ -4.49sro.313. 1.89le~°-° 450s 

6 v ' ' (0.0630) ( 1 .09) [0.248, 1.58] 

3 . 


(5.3) 


8 cg , . _ -0, 0510(0. 086OHO. 818) (48. 0)e 
6 {S) = (0.0630) ( 1 .09) [0.248, 1.58] 


-0.0260s 


(5.4) 


The lateral/directional transfer function fits are plotted as 
dashed lines for comparison with the flight data of Figs. 5.6 and 5.7. 
The matching of the roll-response magnitude (Fig. 5.6a) is excellent 
over the entire frequency range. Only a slight anomaly in the phase 
matching (Fig. 5.6b) is apparent for low-frequency inputs. The agree- 
ment between the sideslip fit and the flight data is also excellent 
(Fig. 5.7). 

As a result of the simultaneous fitting of the roll and yaw 
responses, the transfer-function parameters of Eqs. (5.3) and (5.4) have 
clearly retained their physical significance. The inverse roll mode 
time-constant 1/T r = 1.09 rad/sec is roughly equal to the frequency 
for -45 deg of phase lag as expected. The equivalent Dutch-roll mode is 
lightly damped, with a natural frequency roughly corresponding to the 
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Fig. 5.6. Transfer-Function Model Identification for Roll-Rate 
Response to Aileron (p/-6 a ). (a) Magnitude; (b) Phase. 

peak in the sideslip Bode magnitude plot (Fig. 5.7a). The small effec- 
tive time delays for the roll and sideslip responses reflect negligible 
high-frequency dynamics, and supports the omission of corrections for 
sideslip sensor dynamics. 

The ratio of the natural frequency of the numerator complex zero, 
to that of the denominator complex pole, u> d , is one measure of 
roll-yaw coupling in response to aileron inputs. When the 

numerator and denominator quadratic factors roughly cancel, and the 
resulting decoupled roll response is characterized entirely by the roll- 
mode time-constant. This case typically leads to the best handling 
qualities for a nominal value of the roll-mode time-constant. As the 
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Fig. 5.7. Transfer-Function Model Identification for Sideslip 

Response to Rudder (6 /-<5 ). (a) Magnitude; (b) Phase. 

eg r 

roll-yaw coupling increases, the numerator and denominator quadratic 
factors of Eq. (5.1) no longer cancel, and an undesirable oscillatory 
component of roll rate is generated. With reference to Eq. (5.3), the 
near-unity value of = 1.19 suggests no such concern for roll-yaw 

coupling. Therefore, a simple, decoupled first-order roll-rate model 
could be adopted for future studies. 

The poles and zeros of the identified roll-response transfer func- 
tion are shown on the complex plane in Fig. 5.8. Roll-response handling 
qualities are specified in terms of the location of the dominant pole 
(1/T ), with the Level I (desirable handling qualities) boundary 
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Fig. 5.8. Location of Roll-Rate Transfer-Function (p/6 ) 

3 . 

Parameters in the s-Plane. 

(Ref. 17) shown in the figure. These results indicate very good lateral 
handling qualities for the cruise flight, in contrast to the poor open- 
loop lateral handling qualities for hover. 

The identified transfer functions of Eqs. (5.3) and (5.4) are 
compared with the simulation perturbation transfer functions in 
Table 5.1. Comparison is seen to be excellent for all of the transfer 
function parameters, including the high-frequency sideslip response 
zeros. This comparison is a reflection of the very close frequency- 
response results presented in Figs. 5.2 and 5.4. The slightly more 
shallow phase roll-off for both the roll and sideslip responses 
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(Figs. 5.6b and 5.7b) is reflected in the Dutch roll damping ratio 
discrepancy (11%). 

5.2 Longitudinal Dynamics 

5.2.1 Frequency-Response Identification 

The elevator surface to pitch rate open-loop frequency response 
q/6 0 is shown in Fig. 5.9. The short-period-mode excitation causes a 
peak in the pitch-rate response at about 2.0 rad/sec, with an associated 
phase lag of -45 deg. At higher frequencies, the magnitude response 
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Fig. 5.9. Pitch-Rate Response to Elevator (q/-6 e ). (a) Magnitude; 

(b) Phase. 
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rolls off at -20 dB/decade and the phase shift approaches -90 deg owing 
to the K/s pitch-rate characteristic and negligible flexibility/ 
servolag effects. The drop in the magnitude response and associated 
positive phase response for frequencies below 0.3 rad/sec are due to the 
phugoid dynamics. 

The coherence function shown in Fig. 5.10 is strong over the 
frequency range of 0. 2-7.0 rad/sec. For input frequencies above 



2 

Fig. 5.10. Coherence Function (y r ) for Pitch-Rate Response 

o q 
e M 

Identification. 

7.0 rad/sec, the coherence function becomes erratic and the transfer- 
function identification is less accurate. For low-frequency inputs 
(less than 0.3 rad/sec), the pitch-rate response decreases, even for the 
nearly constant input amplitude, owing to the effect of the phugoid 
dynamics. This results in a decrease of information transfer and an 
associated drop in the coherence function. This coherence function 
roll-off is also attributable to atmosphere turbulence effects, which 
become more important at low frequency where the turbulence spectrum 
peaks (Ref. 46). 
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The pitch-rate frequency response obtained from the simulation 
perturbation transfer function is shown with the flight data in 
Fig. 5.9. The comparison is excellent in the frequency range of 
0. 3-1.0 rad/sec showing that the mathematical model correctly predicts 
the steady-state pitch-rate response to elevator inputs. The differ- 
ences at very low frequency are not conclusive because of the poor 
coherence for frequencies below 0.2 rad/sec (Fig. 5.10) and the lack of 
governor dynamics in the simulation perturbation transfer-function 
models. In the frequency range of 1-10 rad/sec there is a significant 
shift in the magnitude and phase curves between the aircraft and simula- 
tion responses. This discrepancy indicates that both the high-frequency 
gain and dominant response frequency are in error in the model. How- 
ever, the fact that the steady-state gain is correct suggests a dynamic 
response rather than a static response error. This observation in 
addition to the fact that very good static trim correlation has been 
achieved in cruise (Ref. 5), suggests an error in the pitching moment- 
of-inertia, as was also indicated in the hover results. As previously 
discussed, ballast which is used in the aircraft for c.g. control, is 
not included in the present simulator data base and is a likely cause of 
this moment-of- inertia error. The nearly parallel shift between the 
simulation and flight data curves indicates that except for this error, 
the dominant vehicle dynamics are accurately modeled. 

The normal-acceleration response to elevator, a z /6 e , is shown in 
Fig. 5.11. The response is dominated by the classical second-order 
short-period mode over most of the frequency range. The magnitude curve 
is flat at mid-frequency, indicating a constant normal acceleration 
(load factor) response to a step elevator input, with a roll-off in 
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Fig. 5.11. Vertical-Acceleration Response to Elevator (a /6 ). 

(a) Magnitude; (b) Phase. 

response for frequencies beyond the short-period mode. The fall-off in 
normal acceleration response for frequencies below 0.3 rad/sec is prob- 
ably a result of the dominance of the phugoid dynamics. The phase curve 
exhibits the classical second-order characteristic, that is, 0 deg of 
phase lag at low frequency, 180 deg of phase lag at high frequency, and 
90 deg of phase lag at the second-order mode (w = 2.0 rad/sec). This is 
consistent with the previous pitch-rate results. The coherence function 
for the normal acceleration response shown in Fig. 5.12 is strong over 
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Fig. 5.12. Coherence Function (y ) for Vertical-Acceleration 

e z 

Response Identification . 

the wide frequency range of 0.3-10.0 rad/sec, with the fall-off at low 
frequency again resulting from the dominance of phugoid dynamics and 
turbulence effects. The simulation frequency-response is shown with the 
flight data in Fig. 5.11. As in the pitch case, the steady-state 
response to elevator inputs is accurately modeled in the simulation 
(0. 3-1.0 rad/sec). Also, there is a parallel shift at high frequency in 
both the magnitude and phase curves, which as with the previous pitch- 
rate results, indicates an error in the assumed pitching 
moment-of- inertia. 

5.2.2 Transfer-Function Models 

Examination of the longitudinal frequency responses of Figs. 5.9 
and 5.11 indicates that the longitudinal dynamics for this flight condi- 
tion are dominated by short-period mode (w = 2 rad/sec). Therefore, 
the classical pitch-rate and center-of-rotation normal-acceleration 
responses to elevator are applicable: 
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(5.5) 
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where the q(s)/6 e (s) is the Laplace-transformed pitch rate response to 

elevator surface deflection (deg/sec/deg-elevator ) ; M g is the elevator 

pitch sensitivity; 1/T 0 is the first-order numerator inverse time- 
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constant; c gp and u) gp are the equivalent short-period mode damping and 
natural frequency, respectively; and t q is the effective time delay. 
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where a z (s)/6 e (s) is the Laplace-transformed vertical (positive down- 


ward) acceleration response at the instantaneous center-of-rotation to 

elevator surface deflection (g/deg-elevator ) , a„ is the elevator 

% 

vertical sensitivity and the denominator parameters are identical to 
those of Eq. (5.5). The effective time delay for the vertical- 
acceleration response is x a . The low-frequency phugoid dynamics are 
ignored, since, as seen in Figs. 5.9 and 5.11, these are important only 
for inputs of the very lowest frequencies. 

The transfer-function parameters of Eqs. (5.5) and (5.6) are itera- 
tively varied in order to obtain a best least-squares fit between the 
equations and the frequency-responses of Figs. 5.9 and 5.11 over the 
selected frequency range; commonality of the denominator parameters is 
imposed. For a first-cut model, the center-of-gravity (c.g.) and 
longitudinal instantaneous center-of-rotation (ICR) are assumed to be 
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coincident. The iterative fitting procedure is completed using the 
LONFIT computer program. For the present case, the selected range for 
simultaneous fitting of the pitch-rate and normal acceleration responses 
is 0.3-7 rad/sec. In this frequency range, the dynamics are clearly 
dominated by the short-period mode and the coherence is strong for both 
transfer functions. Once the short-period damping and natural fre- 
quency r and a) are obtained for the simultaneous fit, the high- 

frequency gain for the normal acceleration response a is varied, 

5 e 

holding the damping and frequency constant. This procedure is analogous 
to that used in the lateral/directional case and optimizes the a z fit 
over the frequency range 0.3-10.0 rad/sec, since the coherence function 
for this measurement remains strong out to higher frequencies. The 
transfer-function parameters for the pitch-rate and normal-acceleration 
responses are finally obtained: 



-7.73(1 .04)e 


-0.0160s 


(0.554, 2.18] 


(5.7) 


a * cn -0.01 80s 

z , . _ 1 ,60e 

& q (s) = [0.554, 2.18] 


(5.8) 


A comparison of these lower-order models with the flight-test results 
is presented in Figs. 5.13 and 5.14. The pitch-response model 
[Eq. (5.7)] matches the data very well over the selected fitting range 
(0. 3-7.0 rad/sec for q/6 e ), which shows that the short-period approxi- 
mation adequately represents the high- and mid-frequency dynamics. For 
low-frequency inputs (a> < 0.3 rad/sec), the phugoid mode causes a drop 
in the pitch response, a characteristic which is not "captured" by the 
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FLIGHT DATA 

MODEL WITHOUT ICR CORRECTION (EQN. 5.7) 

MODEL WITH ICR CORRECTION (EQN. 5.10) 




Fig. 5.13. Transfer-Function Model Identification for Pitch-Rate 
Response to Elevator (q/-6 e ). (a) Magnitude; (b) Phase. 
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FLIGHT DATA 

MODEL WITHOUT ICR CORRECTION (EQN. 5.8) 

MODEL WITH ICR CORRECTION (EQN. 5.12) 



Fig. 5.14. Transfer-Function Model Identification for Vertical 
Acceleration Response to Elevator (a z /« e ). (a) Magnitude; (b) Phase. 


short-period model. The match between the normal acceleration transfer 
function [Eq. (5.8)] and the flight data is not nearly as good in the 
fitting range, and shows a noticeable discrepancy in the magnitude 
response for frequencies greater than 2 rad/sec. 

Such an error in the mid- and high-frequency a z model could arise 
if the c.g. and ICR are not coincident. Then the vertical acceleration 
response at the c.g. is given approximately by 


-^( 3 ) ; 

e e 


x og s ?;<»> 


(5.9) 
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where x C g is the (unknown) axial location of the center-of-gravity 
relative to the instantaneous center-of-rotation, 

a z 

lCR 

— - — (s) = acceleration response at the ICR 
e 

a 

= ~(s) of Eq. (5.6) 

°e 

3-(s) = Eq. (5.5) 


The frequency responses of Figs. 5.13 and 5.14 were refitted with the 
models of Eqs. (5.5) and (5.9), with \ included in the parameter 
search. The resulting model of pitch-rate response [Eq. (5.5)] is 


_ n ?a(r\ aoro« -0 - 0 °500s 

. -7.3o(0.o90)e 

6 e ’ [0-536, 2.02] 


( 5 . 10 ) 


The vertical acceleration response at the center-of-rotation [Eq. (5.6)] 
is 



1.17 

[0.536, 2.02] 
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(5.11) 


and the location of the c.g. relative to the ICR is x na = -5.73 ft 

eg 

(c.g. behind ICR). The resulting acceleration response at the c.g. is 
obtained from Eq. (5.9): 


_fc g, ) _ -0 . 0230( -6 . 70 ) ( 7 . 59 ) e'° ‘ 00200s 
6 e [0.536, 2.02] 


(5.12) 
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which is a higher numerator-order model than the simple model of 
Eq. (5.8), and displays the nonminimum phase (right-half plane) zero 
that results from the c.g. being located behind the ICR. 

The pitch-rate model of Eq. (5.10) is shown in Fig. 5.13 for com- 
parison with the previous results. Neither the frequency response nor 
the transfer-function parameters of the q/6 0 model are significantly 

altered by including the x n <r degree-of-freedom. However, the c.g. 

c a 

acceleration model of Eq. (5.12) fits the extracted frequency response 
much more closely, as shown in Fig. 5.14. Thus the noncoincident c.g. 
and ICR appear to be a likely source of the discrepancy obtained when 
the simple model of Eq. (5.8) is used. 

The identified short-period damping and frequency (5 ar> = 0.536 and 

s>p 

u> = 2.02 rad/sec) are consistent with preliminary observations, which 
sp 

were based on the raw Bode plot information. The relatively high degree 
of damping is a reflection of the small response peak of the pitch-rate 
transfer function. The small effective time delays (in the pitch atti- 
tude and normal acceleration transfer-function fits) indicate that the 
high-frequency flexibility effects are not important for this flight 
condition, as noted earlier. The military standards requirements 
(Ref. 17) for lower-order equivalent system pitch response are given in 
terms of the parameters w S p T e 2 > t S p> t q . The values of these parame- 
ters given in Eq. (5.10) are well within the Level I, category A, indi- 
cating desirable longitudinal handling for cruise flight. 

The identified transfer functions are compared with the simulation 
perturbation transfer functions in Table 5.1. The simulation transfer 
functions contain an additional complex mode associated with the phugoid 
dynamics and a nearly canceling pole/zero pair associated with the rotor 
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angular velocity (rpm) degree-of-freedom. The fact that the rpm mode 
nearly cancels and that the phugoid mode is at very low frequency justi- 
fies ignoring these dynamics and adopting the simple short-period 
models. 

The pole and zero locations of the simulation transfer functions 
compare very well with the flight-extracted results. The good corre- 
spondence of the numerator factors in the vertical acceleration response 
supports the corrections based on instantaneous center-of-rotation 
effects [Eqs. (5.9), (5.12)]. The differences in short-period frequency 
and pitch-response sensitivity can both be explained by the error in 
pitching moment-of- inertia. The pitch-axis sensitivity is the 

pitching moment (foot-pounds) per inch of elevator normalized by the 
pitching inertia (Iy) (Ref. 41). The ratio of pitch sensitivities for 
the simulation and flight results (Table 5.1) gives: 


(M 6 > 
e sim 

(M > 

e flight 


= 1.67 = 
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The short-period frequency ui depends on the square-root of 

sp 

stability derivatives which are also normalized by the pitching moment- 
of-inertia (Ref. 41). Based on the short-period frequency: 


t*— = 1.29 = Vi. 67 = UfH 

w sp flight L sim J 


(5.14) 


which is the same result as in Eq. (5.13). Thus, an error in pitch 
inertia would explain the differences in both pitch sensitivity and 
short-period frequency, while not disturbing the good static trim cor- 
relation. An underestimation of pitch inertia in the simulation is 
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indicated from both the hover and cruise results, although there is some 
discrepancy in the error ratio (67? from the cruise results and 25? from 
the hover results). 

5.3 Verification of Transfer-Function Models 

The high degree of inherent bare-airframe stability for the cruise 
flight condition makes open-loop verification inputs much easier to 
complete than in the hover condition. Steady initial conditions and 
clean long-duration (6-8 sec) step inputs are easily obtained. The 
aircraft (filtered) roll-rate response to a step aileron input is shown 
in Fig. 5.15. The transfer- function model of Eq. (5.3) shown in the 

FLIGHT DATA 

MODEL 




Fig. 5.15. Roll-Response Model Verification Using SCAS-off Flight 
Data, (a) Aileron Input; (b) Roll-Rate. 
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dashed line matches the flight data quite well for the entire input 
duration, except for some discrepancy in steady-state roll rate. This 
error may be associated with the p/6 a transfer-function model mismatch 
at mid-frequency (Fig. 5.6), some nonlinear effects associated with the 
large steady-state roll rate, or simply a gust. The sideslip response 
to a step rudder input is shown in Fig. 5.16, with the appropriate 
corrections for the nose-boom position. The transfer-function model of 
Eq. (5.4) matches the flight data almost identically for the entire 
time-history. 


FLIGHT DATA 



0 2.5 5.0 7.5 10.0 12.5 

TIME, sec 

Fig. 5.16. Sideslip-Response Verification Using SCAS-off Flight 
Data, (a) Rudder Input; (b) Sideslip Response. 


For most of the flight tests, frequency-sweeps and step inputs were 
performed consecutively for each flight condition, thus maintaining 
constant altitude, airspeed, and loading condition. Unfortunately, the 
open-loop elevator step response data for = 170 knots were taken 
at H d = 2800 ft (density altitude), whereas H d = 10,000 ft for the 
sweeps. Since the indicated airspeed was constant and, therefore, also 
the dynamic pressure, the derived pitch-rate model is valid for the 
lower-altitude step responses. However, the significant difference in 
true airspeed (about 20 knots) and associated centripetal acceleration 
(= U Q q) makes an a z step response comparison invalid for these data. 
The response to elevator doublet inputs, obtained from closed-loop 
(SCAS-on) tests, for = 170 knots and = 12,000 ft is used to 
validate the a z model, although the aircraft loading condition (and 
altitude) is still not exactly matched. 

The SCAS-off (open-loop) pitch-rate response to step elevator 
deflection (H^ = 2800 ft) is shown in Fig. 5.17. The elevator input 
corresponds to about *10% of the maximum control deflection. The compar- 
ison between the ICR-corrected model response [Eq. (5.10)] and the 
flight data is seen to be very good over most of the time-history. The 
slight deviations occurring toward the end of the run are due to the 
inadequacy of the short-period approximation in modeling the low- 
frequency (phugoid) dynamics (see Fig. 5.13a). Even so, the short- 
period approximation clearly gives an excellent characterization of the 
important initial response dynamics. 
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Fig. 5.17. Pitch-Response Model Verification Using SCAS-off Flight 
Data (h d = 2800 ft), (a) Elevator Input; (b) Pitch Rate. 

The pitch-rate model also compares favorably with the flight data 
for the response to an elevator doublet input, as shown in Fig. 5.18c. 
However, the normal-acceleration model of Eq. (5.12) is not very good. 
The initial reversal predicted by the model (effect of the x C g correc- 
tion) does not occur in the flight response in Fig. 5.18c. In fact, if 
the original center-of-rotation model of Eq. (5.8) is used, the match 
with the flight data is much better, as is also shown in Fig. 5.18c. 

These results indicate that the a / 6 model which includes the 

z cg e 

center-of-rotation correction [Eq. (5.12)] is highly tuned to the 
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measured frequency response for sinusoidal inputs (Fig. 5.14) and is not 
robust to small changes in flight condition and input form. The fact 
that acceleration is a higher-order state derivative than the other 
angular-rate derivatives further exaggerates the sensitivity of this 
model, especially for the initial (high-frequency) response dynamics. 
Therefore, the lower-order model of Eq. (5.8) should be used for future 
handling qualities and control-system analyses of this flight condi- 
tion. This is a good example of the dangers of model overparameteriza- 
tion. The capability to evaluate directly the effect of model order on 
response fitting is also clearly seen in this example, again represent- 
ing a unique advantage of the frequency-domain approach for this kind of 
comparison. 

Overall, these verification results indicate that the derived 
lower-order transfer functions are very good predictors of the open-loop 
vehicle response in cruise over a wide frequency range. 
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Chapter 6 
CONCLUSIONS 


The following conclusions are separated into those specifically 
associated with the analysis of the XV- 15 and those associated with the 
analysis of the frequency-response identification method in general. 

6.1 Analysis of the XV- 15 

Specific findings from the identification of the open-loop hover 
dynamics are as follows. 

1 . Heave response to power lever and yaw response to rudder are 
decoupled from the other degrees-of-freedom and are very lightly damped. 

2. Pitch and roll open-loop dynamics exhibit unstable low- 
frequency oscillations. High pilot work-load should be expected in the 
event of a SCAS failure in either of these axes. 

3. Significant coupling of yaw response to aileron has been iden- 
tified. Pedal deflection equal to about 50 % of lateral stick inputs is 
necessary to maintain a fixed heading. 

4. The effective high-frequency time-delays are small in all axes 
and are probably not important in determining the handling-qualities of 
the XV- 15 aircraft. 

5. The comparison of piloted simulation and flight-test dynamic 
results for high-frequency inputs (w > 1.0 rad/sec) is generally favor- 
able in all axes, with larger anomalies existing for low-frequency 
inputs. However, underestimates in the pitch and yaw inertias used in 
the simulation model caused the overprediction of pitch- and yaw-control 
sensitivities. 
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6. Small-perturbation transfer functions obtained from the 
piloted-simulator (VMS) are poor representations of the low-frequency 
dynamics. However, perturbation transfer functions for the current 
nonreal-time simulation compare more favorably with the flight data, 
suggesting that the corrections to the mathematical models have been 
appropriate. 

7. The step responses of the identified transfer-functions and the 
aircraft generally match quite well, even for large-amplitude motions. 
This verifies that an accurate representation of the open-loop hover 
dynamics has been achieved, and that the hover transfer-function models 
derived herein can be used for future XV-15 control system and handling- 
qualities analyses. 

Specific findings from the identification of the open-loop cruise 
(V^ = 170 knots) dynamics are as follows. 

1. The response characteristics are very stable and decoupled, and 
indicate desirable open-loop handling-qualities. 

2. Standard lower-order equivalent models adequately match the 
identified dynamics in all degrees-of-freedom except for the normal 
acceleration responses to elevator deflection. Although a center-of- 
rotation correction is needed to accurately fit the a z / fi e frequency- 
response data, the resulting higher-order model is not robust to small 
changes in flight condition and input form. 

3. Small perturbation transfer functions and associated frequency 
responses from the nonreal-time simulation correlate very well with the 
flight-extracted results except for an error in the simulation pitch and 
yaw inertias. 
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4. The step responses of the identified transfer functions and the 
aircraft generally match very well. 

6.2 Frequency-Response Identification Method 

1 . Non-parametric frequency-response identification techniques 
produce an excellent physical understanding of basic vehicle dynamics. 

2. Non-parametric techniques are especially useful for documenting 
the response characteristics of a new configuration for which a prior 
parametric model is not known. 

3. Extraction of open-loop responses from closed-loop data is 
feasible with frequency-sweep inputs. Closed-loop testing of rotorcraft 
with highly unstable bare-airframe dynamics allows the execution of long 
duration inputs needed for low-frequency response identification. 

4. Multi-input/multi-output frequency-domain methods are suitable 
for extracting the coupled responses of rotorcraft. 
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Chapter 7 

RECOMMENDATIONS FOR FUTURE STUDY 


The objectives of the XV-15 study (Sec. 1.1) have been achieved. 

The following recommendations pertain to future studies of the 
frequency-response identification method. 

Future studies in frequency-response identification should focus on 
advancing the method in two major areas: 

1. Reconstruction of physical stability-derivatives from the 
transfer-function matrix 

2. Generalization of two-input/single-output frequency-response 
identification to multi-input/multi-output (MIMO) systems 

7.1 Reconstruction of Physical Stability Derivative Parameters 

The frequency-response method illustrated in this study yields a 
transfer-function matrix for the coupled system. The next logical step 
in a full identification procedure is to reconstruct the physical 
stability-derivative (state-space) realization rather than a canonical 
realization as in the present aproach. The MIMO state-space reconstruc- 
tion approach given by Chen (Ref. 47) is well suited for this purpose 
because it takes into account the physical redundancy in the transfer- 
function zeros. Frequency-domain parameter identification would thus be 
completed in a three step procedure: 

Step 1 . Frequency-response identification 

Step 2. Transfer-function modeling 

Step 3. Stability-derivative reconstruction 
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By breaking the parameter-identification approach into these three 
distinct steps, maximum benefit can be gained from the nonparametric 
identification before the parametric step is attempted. This approach 
may yield both a more accurate, and physically meaningful and consistent 
state-space model than can be identified using direct time-domain or 
frequency-domain parameter-identification procedures. 

Stability-derivative reconstruction should be tested on the XV- 15 
transfer-functions for the hover flight condition, since a complete 
state-space model has previously been identified for this condition 
using time-domain identification methods (Ref. 22). The frequency- 
domain parameter-identification procedure should then be tested on a 
single-rotor helicopter which has more highly-coupled system dynamics 
than the XV- 15. 

7.2 Generalization to MIMO Systems 

Spectral quality assessment is conducted on each input/output 
record pair to determine the runs suitable for concatenation and best 
window sizes. Although this assessment is relatively straightforward, 
the process is computer and labor intensive, and could be automated 
using a rule-based algorithm. This would be especially useful for MIMO 
identification where the number of required frequency responses 
increases rapidly with additional cross-coupling paths. 

The two-input identification software tools developed (or adapted) 
in this study (TISOSA, TFTHISTORY, NAVFIT, LONFIT, LATFIT) need to be 
further improved to allow MIMO identification and verification of sys- 
tems with more than two control inputs. A general matrix formulation 
will allow extension to the multi-input case. 
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Appendix A 
CHIRP z-TRANSFORM 


Existing computer programs for identifying frequency-responses from 
flight data generally use the standard (Cooley-Tukey) FFT procedure to 
transform data from the time-domain to the frequency-domain (Refs. 10, 

11, 13). Initial versions of the FRESPID program also used a standard 
FFT algorithm, but this was later replaced by the chirp z-transform 
(CZT). The mathematical basis for this algorithm is given in Refs. 28 
and 32, and is not discussed in this appendix. The CZT provides more 
flexibility in selection of computational parameters than is allowed by 
the standard FFT. However, it has increased computational time and 
storage requirements. 

Three important differences between the FFT and the CZT are 
reviewed in this appendix: (1) Resolution, Af; (2) Starting frequency, 

f0; and (3) Window size, LAt. These parameters define the discrete 
Fourier transform calculation ( DFT ) , which is repeated from Eq. (3.2): 


N-1 

X(f. ) = X(kAf) = At x exp[-j2w(kn)/N] ; k = 0, 1, 2,..., 
K n=0 n 

where 

X(f k ) = Fourier coefficients 
x n = x(nAt) 

At = time increment 

N = number of discrete frequency points 


N-1 


( A. 1 ) 
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A.1 Resolution 


The (frequency) resolution of the standard FFT (Af) is determined 
by the width of the window (LAt, sec), since for this algorithm the 
number of time-history points (L) is also the same as a number of output 
spectral points (N): 


Af = 


LAT 


Hz 


(A. 2) 


which in terms of the sampling frequency f g = 1/At is 

f f 

Af = ^ ^ (A. 3) 

Thus, the FFT output frequencies are N evenly spaced points 

around the unit circle. Since the z-transform is a mirror image with 

respect to the Nyquist frequency (f = f /2), the second half of the 

c s 

spectrum (f /2 < kAf < f ) is redundant and is not actually calculated 
s s 

in the FFT procedure of Eq. (A.1). Furthermore, the Nyquist frequency 
is generally selected to be much greater than the upper end of the 
identification frequency range (0.1-10.0 rad/sec = 0.016-1.59 Hz in the 
present study) to avoid digital distortion, so many of the spectral 
calculation points are not needed (1.59 < kAf < f /2). 

The basic difference in the CZT algorithm is that the spectral 
calculation is completed along an arbitrary specified arc in the 
z-plane. When, as in the present case, this arc is taken as a segment 
of the unit circle, the CZT provides the Fourier coefficients. A key 
characteristic of the CZT is that the number of calculation points N 
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within this segment is arbitrary; thus, the frequency resolution is also 
arbitrary'!’ : 



(A. 4) 


where 

f Q = starting spectral frequency (arbitrary, Sec. A. 2) 
f.j = ending spectral frequency (arbitrary, Sec. A. 2) 

A commonly held misconception is that arbitrary frequency resolu- 
tion can also be achieved using the standard FFT by simply "padding" the 
input sequence with zeroes (artificially increasing L). However, pad- 
ding with zeroes only allows interpolation between true frequency reso- 
lution points (Af = f g /L) and does not increase the frequency resolu- 
tion . The only way to increase the true frequency resolution of the 
standard FFT is to increase the width of the window, which decreases the 
amount of averaging and so increases the random error (Sec. 3.2.3). 


A. 2 Starting Frequency 

The minimum output frequency f Q that can be identified using the 
standard FFT procedure is also the frequency spacing: 



(A. 5) 


t The author has encountered numerical problems with the CZT algorithm 
when the number of output points (N) is much greater than the number of input 
points (L), N 2 2L. However, since all the output points are concentrated 
within the desired frequency range, the resulting frequency resolution is 
still substantially higher than is achievable with the standard FFT procedure. 
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Reducing the starting frequency to allow an identification of the 
lower frequency dynamics means selecting wider windows and, therefore, 
reducing the amount of spectral averaging [and increasing the spectral 
variance, Eq. (3.10)]. In the CZT, the initial (and final) frequency 
f Q (and f.,) is arbitrary. Thus, the initial frequency is independently 
selected to correspond to the lower end of the frequency range of inter- 
est (fQ = 0. 1 /2n = 0.016 Hz in this case). 

ft. 3 Window Size 

In both the standard FFT and CZT procedures, window size (LAt) 
controls the balance between low-frequency signal content and spectral 
averaging (Sec. 3.2.1). In the standard FFT, the number of time-history 
points in the window must be an integer power of 2 (a "composite 
number," e.g., 512, 1024, 2048). Thus, the allowable adjustment in the 
width of window is coarse (only halving or doubling in size), for a 
fixed sample interval (At). The chirp z-transform allows the arbitrary 
selection of L (including prime numbers). The capability to indepen- 
dently set the starting frequency (f Q , Sec. A. 2) and then finely tune 
the window size yields both improved coherence and lower spectral vari- 
ance over a wide frequency range as compared to the standard FFT . 

These three characteristics of the chirp z-transform make it a 
significantly more flexible algorithm than the standard FFT; this flexi- 
bility can yield a substantial improvement in the spectral resolution 
and dynamic range, especially at the low-frequency end. The author has 
found the flexibility of the chirp z-transform to be very important in 
achieving high-quality frequency-response identification from (noisy) 
flight data. 
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Appendix B 

LIMITATIONS IN CLOSED-LOOP IDENTIFICATION 


The identification of low-frequency helicopter dynamics requires 
long-duration inputs. In hover flight the dominant open-loop dynamics 
are unstable, so execution of these inputs with the stability and con- 
trol augmentation system disengaged (SCAS-off) is not practical. There- 
fore, it is desirable to extract the open-loop dynamics from flight-test 
data obtained with the augmentation system engaged ("closed-loop identi- 
fication"). 

This appendix derives the fundamental relationships for single- 
input/single-output closed-loop identification in the presence of pro- 
cess noise. A simple relationship relating bias errors in the identi- 
fied open-loop frequency-response to the input-to-process noise ratio is 
developed. A numerical study illustrates the important concepts using a 
simple model of the XV- 15 unstable dynamics in hover. 

B.1 Spectral Relationships 

Consider the identification of the single- input/single-output roll- 
rate response to aileron from closed-loop tests as illustrated in 
Fig. B.1. The closed-loop system is the same as that previously dis- 
cussed in Sec. 3.2.4 (Fig. 3.7), except in the present case the compen- 
sation element (G c ) is placed in the forward loop. This provides a 
steady-state closed-loop gain of approximately unity — independent of the 
compensation gain. The other elements of Fig. B.1 are as follows. 

1 . 6^ = pilot lateral-stick input 

2. p = measured roll-rate (measurement noise is neglected here) 
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Fig. B. 1 . Closed-Loop Identification. 


3. 6 = total aileron surface deflection 

4. H = pg/6 & = open-loop roll-rate frequency-response 

5. n = roll-rate process noise. 

Some common sources of noise are: (a) turbulence, the auto- 

spectrum of n is determined from the turbulence velocity autospectrum 
and the roll-rate response to turbulence velocity inputs; (b) neglected 
secondary inputs; (c) nonlinear remnants which are not accounted for by 
the linear open-loop frequency-response H. The external inputs from 
the pilot and process noise sources are assumed to be uncorrelated; 

viz. , G„ =0. 

6 L n 

The desired open-loop frequency-response estimate is obtained from 
Eq. (3.8): 


The cross-spectrum Gg p is determined from Fig. B.1 (Ref. 28): 


G v = °v. + V 
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Substituting Eq. (B.2) into Eq. (B.1): 


H = 


6 a P * 


6 n 
a 


6 6 
a a 


(B.3) 


H 


6 n 
a 


6 6 
a a 


True response Bias error 


(B.4) 


which shows that the spectral estimate is biased by correlation between 
the aileron surface and process noise inputs. This bias term is zero 
for open-loop systems since the pilot and process noise inputs are 
uncorrelated; the undesirable correlation results from the feedback of 
the process noise around the loop. Equation (B.4) also shows that the 
true (unbiased) open-loop frequency-response is identified regardless of 
the level of compensation when process noise is not present. 

The cross-spectrum is determined from the closed-loop rela- 

tionship between process noise inputs and aileron surface response 
(Ref. 28)1": 


* 



1" Asterisk (*) denotes complex-conjugate. 
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The aileron surface input-autospectrum is composed of components 


from the two external (uncorrelated) excitations: 

2 


G S 6 G 6 r 6r 
3 , 3 , L Li 


6 

2 

6 

a 


a 


+ G 


6 j 

nn 

n 

L 




= G 


6 l 6 l 


G 

2 

-G 

c 

, n 

c 

1 + G H 

+ U 

nn 

1 + G H 

c 


c 


Substituting Eqs. (B.5) and (B.6) into Eq. (B.4): 


H = H - 


G [G /(I + G H)] 
nn c c 


|G /(I + G H) | \G + G ) 
1 c c 1 nn 


But, 


V • (t \) (A?) 


SO, 


H = H 


nn 


t G 0 /(1 * 0 o H)][G * G nn ] 

L L 


( B . ( 


(B. i 


(BA 


(BA 


When the closed-loop system is excited only by process noise (i.e 


no pilot inputs), G =0 and 

°l°l 


nn 


[G /(I + G H) ]G 


H = H 



which is the limiting condition of identifying the inverse compensa- 
tion. This result is in agreement with the simple analysis of 
Eqs. (3.11 and 3.12). 

When pilot and process noise inputs are both present , we can define 
the signal-to-noise ratio: 


2 \\ M 

s M s c <«,) 

nn 


(B. 1 1 ) 


which depends on frequency. 

Substituting Eq. (B.11) into Eq. (B.9) and normalizing the estimate 
by the true frequency-response yields the following normalized bias 
error: 


. b[»l = zl 

b H [G H/(1 + G H)][S 2 (w) + 1] 


(B. 12) 


where it should be recognized that the denominator term in Eq. (B.12) 
G C H/[1 + G C H] is the closed-loop transfer-function. The open-loop 
unstable modes of rotorcraft are generally at frequencies well below the 
open-loop crossover frequency u> c (defined as frequency at which 
|G C H| = 1) for hovering vehicles. In the case of the XV-15, the domi- 
nant open-loop unstable modes are at a frequency of about 
to = 0.5 rad/sec [Eq. (4.4)], compared to a crossover frequency of 
u> c = 2.5 rad/sec (closed-loop time-constant = 0.4 sec, Sec. 2.1.4). At 
low frequency (w < a> c ), the closed-loop transfer-function becomes 
(Ref. 41) 
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1 


(B. 13) 


G H 
c 

1 + GH 
c 


0)< CD 


c 


which when substituted into Eq. (B.12) gives the following simple 
expression for the normalized bias error in the frequency range of the 
open-loop dynamics: 


_ 1 
b S 2 (w) + 1 


(B. 14) 


Equation (B.14) indicates that the normalized bias error can be 
kept to within 10% (e^ < 0.10) for a signal-to-noise ratio: 

S(u») > 3 (B. 15) 


In terms of the noise-to-signal ratio o(u>) = 1/S (w), the requirement is 

8(u>) < 0.33 (B. 16) 

which was quoted in Sec. 3.2.4. 

This analysis shows that the pilot's input must exceed the input 
from process noise by at least a factor of 3 across the frequency range 
of interest to reduce the normalized bias error to an acceptable 
value. The use of the frequency-sweep, which generates a persistent and 
relatively large vehicle response over the entire identification fre- 
quency range, ensures that this condition is maintained and that the 
bias errors are small. Multi-step inputs generally do not have uniform 
power distribution. If process noise inputs are large in the frequency 
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range where the multi-step signal content is small, large bias errors 
may result, making such inputs poorly suited for closed-loop 
identification. 

B.2 Numerical Study 

A numerical study was conducted to illustrate the concepts dis- 
cussed in the preceding section and to validate the expression for 
normalized bias error given in Eq. (B.14). A simple model of the XV-15 
dynamics in hover is used to facilitate a clear physical understanding 
of the limitations of closed-loop identification, which is important for 
future applications. 

B.2. 1 Simplified Closed-Loop Model 

The open-loop roll-rate response in hover [Eq.(4.4)] is represented 
by a simple (normalized) first-order unstable transfer-function: 

-0.005s 

H<s)s f-rccr 

where the 5 msec time-delay is included to account for numerical 
integration. 

The compensation transfer-function is assumed to be a pure gain: 

G = K (B.18) 

c 

The closed-loop time-constant of the XV-15 roll-response in hover 
is 0.4 sec (Sec. 2.1.4), which implies a closed-loop first-order roll- 
mode at s = -2.5 rad/sec, yielding a compensation gain of: 
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K = 3.0 


(B. 19) 


The pilot lateral-stick excitation was derived by passing a 
Gaussian random sequence (200-sec duration) through a low-pass filter 
(1.0 rad/sec roll-off frequency). This yields an input-autospectrum 
(G g g ) approximating that of the lateral-stick frequency-sweep of 
Fig. 4.2. For illustrative simplicity, the process noise spectrum was 
also obtained by passing a Gaussian random sequence through the 
1.0 rad/sec low-pass filter. This yields a signal-to-noise ratio S(w) 
which is independent of frequency, and is just equal to the stick input- 
to-noise rms ratio. 

A simulation program was developed to numerically integrate the 
governing equations of Fig. B.1 and provide time-histories for 
frequency-response identification using the FRESPID program. Results 
were obtained as a function of the noise-to-signal ratio a = 1/S for 
the nominal compensation gain of the XV- 15 (K = 3.0) and a gain typical 
of a high-bandwidth control system (K = 6.5 rad/sec). 

B.2.2 Results of Numerical Study 

The open- and closed-loop frequency-responses (p/6 and p/6, , 
respectively) for the nominal gain (K = 3.0) are shown in Figs. B.2 
and B.4, along with the ideal transfer-functions. The associated coher- 
ence functions shown in Figs. B.3 and B.5 indicate nearly perfect iden- 
tification over the entire frequency range, with a slight random error 
for the open-loop response in the frequency range near the dominant 
open-loop mode (w = 0.5 rad/sec). The open-loop frequency-response is 


150 


IDENTIFICATION 



Fig. B.2. Identification of Open-Loop Frequency Response for the 
Nominal-Gain System, K = 3.0. (a) Magnitude; (b) Phase. 



Fig. B.3. Coherence Function for Open-Loop Roll-Rate 
Identification, Nominal -Gain System, K = 3.0. 
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IDENTIFICATION 



.1 1 10 
FREQUENCY, rad/sec 


Fig. B.4. Identification of the Closed-Loop Frequency-Response for 
the Nominal-Gain System, K = 3.0. (a) Magnitude; (b) Phase. 



FREQUENCY, rad/sec 

Fig. B.5. Coherence Function for Closed-Loop Response 
Identification, Nominal-Gain System, K = 3.0. 


also identified with almost no error for the high-gain case (Figs. B.6, 
B.7) ; the higher-bandwidth closed-loop frequency-response matches the 
true response (Figs. B.8, B.9). 
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FREQUENCY, rad/sec 


Fig. B.6. Identification of the Open-Loop Frequency-Response for 
the High-Gain System, K = 6.5. (a) Magnitude; (b) Phase. 



FREQUENCY, rad/sec 

Fig. B.7. Coherence Function for Open-Loop Response 
Identification, High-Gain System, K = 6.5. 
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Fig. B.8. Identification of the Closed-Loop Frequency-Response for 
the High-Gain System, K = 6.5. (a) Magnitude; (b) Phase. 
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Fig. B.9. Coherence Function for Closed-Loop Response 
Identification, High-Gain System, K = 6.5. 


The effect of increasing the noise- to-signal ratio a on the 
open-loop identification for the nominal gain system is shown in 
Figs. B.10 and B. 11. Very little error is apparent for noise-to-signal 
ratios of o < 0.3. At moderate noise-to-signal ratios (o = 0. 3-1.0), 
the identification errors (associated with both bias and random error 
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Fig. B . 1 0 .* Effect of Noise-to-Signal Ratios a = 0, 0.1, 0.3 on 
the Identification of the Open-Loop Frequency-Response for the Nominal- 
Gain System, K = 3.0. (a) Magnitude; (b) Phase; (c) Coherence. 
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Fig. B.11. Effect of Noise- to-Signal Ratios 3=0, 1.0, 10.0 on 
the Identification of the Open-Loop Frequency-Response for the Nominal- 
Gain System, K = 3.0. (a) Magnitude; (b) Phase; (c) Coherence. 
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contributions) are significant and are reflected by the drop in the 
coherence function. Thus, for low and moderate noise levels, the coher 
ence function remains a good indicator of spectral accuracy for closed- 
loop identification. The case of nearly infinite noise-to-signal ratio 
(5 = 10) shows that the open-loop frequency-response identification 
approaches the inverse compensation transfer-function as indicated in 
Eq. (B. 10) : 



= 0.333 = -9.54 dB 



-180 deg 


(B.20) 


and the associated coherence function approaches unity (Fig. B.llc). 

A better quantitative appreciation for the effect of process noise 
on the identified open-loop frequency-response is achieved by looking at 
two measures of error. The first measure is the total rms error . This 
relates the deviation between the estimated open-loop frequency- 
response H and the true open-loop frequency-response H [Eq. ( B . 1 7 ) 3 : 


[ 20 "| 

20 ^ [H dB (u>) * H dB (u)] J 

°e p = [k £ <<“> - W“» 2 ] 


(B.21 ) 


These calculations are performed by the NAVFIT program, which 
selects 20 linearly spaced frequency points across the logarithmic scale 
(Sec. 3.3.2). The total rms magnitude and phase errors are plotted as a 
function of noise-to-signal ratio o in Figs. B.12a and B.12b for the 
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— •— NOMINAL XV-15 GAIN, K = 3 
— o— HIGH GAIN, K = 6.5 




Fig. B.12. Variation in the Total rms Errors with Noise-to-Signal 
Ratio for the Identification of the Open-Loop Frequency-Response. 

(a) Magnitude Error; (b) Phase Error. 

nominal and high-gain cases. Based on an allowable magnitude error of 

10 % (o 0 = 0.92 dB), the results of Fig. B.12a indicate that the noise- 

g 

to-signal ratio should be a < 0.11 for the low-gain case and 
a < 0.20 for the high-gain case. Adopting the relative error weighting 
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of 1 dB magnitude error to 7 deg phase error (Sec. 3.3.2) yields an 

allowable rms phase error o = 6.4 deg (= 7 * 0.92 dB). Based on this 

P 

phase error limit, the results of Fig. B.12b indicate an acceptable 
noise-to-signal ratio of a < 0.22 for the low-gain case and 
a < 0.25 for the high-gain case. Thus, the noise limits based on the 
magnitude criteria are more restrictive for this study. The larger 
allowable levels of noise for the high-gain case occur because the 
resulting closed-loop system is more resistant to excitation from pro- 
cess noise, thereby reducing the total level of random error. For large 
noise-to-signal ratios o >> 1, the magnitude errors asymtotically 
approach different levels for the two feedback gain cases because the 
associated inverse compensation gain values -1/K are different. The 
rms phase errors approach a common value because the open-loop identifi- 
cation tends toward $ = -180 deg, independent of the gain level 
Eq. (B.20) . 

A direct measure of parameter bias error is obtained by fitting 
each of the frequency-responses in Fig. B.10 (and similarity for the 
high-gain case) with the first-order transfer-function model: 

K p" TS 

H = * 6 o (B.22) 

S “ a 


and tabulating the bias errors normalized with respect to the true open- 
loop transfer-function parameters of Eq. (B.17): 


e b tK] 


e b [a] 


MU 

K 


b[a] 

a 


|K - 1 


a - 0.5 
0.5 


100 $ 


(B.23) 


100 $ 
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The normalized bias errors for transfer-function gain e b [K] and 
mode location e^ta ] are plotted in Figs. B.13a and B . 1 3b . The small 

— •— NOMINAL XV-15 GAIN, K = 3 
— o— HIGH GAIN, K = 6.5 




Fig. B. 13. Variation in Normalized Bias Errors with the Noise-to- 
Signal Ratio for the Identification of the Open-Loop Frequency Response 
(a) Normalized Gain Error; (b) Normalized Mode Location Error. 


bias error for no process noise is associated with the overall spectral 
identification procedure (Ref. 28). The increase in bias error for the 
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high-gain case at low process noise levels (o < 0.1) probably results 
from the integration errors which become more significant at the higher 
closed-loop natural frequencies. Restricting the maximum bias error 
(gain and mode) to within 10# indicates a requirement for noise-to- 
signal ratios of o < 0.3, which is in good agreement with the theoreti- 
cal results of Eq. ( B . 1 6 ) . The fact that the bias errors are very small 
for process noise ratios of o < 0.10 indicates that the random error 
contribution is the significant factor in the total rms error results of 
Fig. B.12a and B.12b. By limiting the noise ratio to o < 0.10 in 
order to maintain acceptable levels of total rms error, negligible bias 
errors are ensured. 

B.3 Conclusions of Closed-Loop Identification Study 

Two major conclusions can be drawn from this study of limitations 
on closed-loop identification: 

1 . Identifying the open-loop frequency-response from closed-loop 
flight data results in a biased estimate whenever process noise is 
present. Common sources of process noise which must be considered are 
turbulence, extraneous unaccounted for inputs, and nonlinearities. 

2. Process noise causes random and bias errors in the estimates. 
Numerical studies based on the simplified model of XV- 15 dynamics 
suggest that allowable levels of noise are determined more from the 
effects of random error than from those of bias error. Total error can 
be kept to an acceptable level when the noise-to-signal rms ratio is 

o < 0.10. Bias errors do not become important until process noise 
ratios reach o > 0.3. The coherence function is a good indicator of 
the level of total random error for low and moderate process noise 
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levels, and can be used to assess the quality of the open-loop 
frequency-response identification from closed-loop test data. 

This relatively simple numerical simulation provides a good appre- 
ciation for the effects of process noise and closed-loop compensation. 
Similar studies should be conducted (perhaps using more sophisticated 
models of vehicle response and noise environment) to support future 
closed-loop identification experiments on rotorcraft. 
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